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ABSTRACT. Further genomic sequencing of butterflies by our research group expanding the coverage of species and 
specimens from different localities, coupled with genome-scale phylogenetic analysis and complemented by phenotypic 
considerations, suggests a number of changes to the names of butterflies, mostly those recorded from the United States and 
Canada. Here, we present evidence to support these changes. The changes are intended to make butterfly classification more 
internally consistent at the genus, subgenus and species levels. I.e., considering all available evidence, we attempt to assign 
similar taxonomic ranks to the clades of comparable genetic differentiation, which on average is correlated with the age of 
phylogenetic groups estimated from trees. For species, we use criteria devised by genomic analysis of the genetic 
differentiation across suture zones and comparison of sympatric populations of closely related species. As a result, we resurrect 
4 genera and 1 subgenus from subgeneric status or synonymy, change the rank of 8 currently used genera to subgenus, 
synonymize 7 genus-group names, summarize evidence to support 19 taxa as species instead of subspecies and 1 taxon as 
subspecies instead of species, along with a number of additional changes. One new genus and one new subspecies are 
described. Namely, the following taxa are treated as genera Tharsalea Scudder, 1876, Helleia Verity, 1943, Apangea Zhdanko, 
1995, and Boldenaria Zhdanko, 1995. Tetracharis Grote, 1898 is a valid subgenus (not a synonym of Anthocharis Boisduval, 
Rambur, [Dumenil] & Graslin, [1833]) that consists of Anthocharis cethura C. Felder & R. Felder, 1865 (Muller, 1764), 
Anthocharis midea (Hiibner, [1809]), and Anthocharis limonea (A. Butler, 1871). The following are subgenera: Speyeria 
Scudder, 1872 of Argymnis Fabricius, 1807; Aglais Dalman, 1816 and Polygonia Hiibner, [1819] of Nymphalis Kluk, 1780; 
Palaeonympha Butler, 1871 of Megisto Hiibner, [1819]; Hyponephele Muschamp, 1915 of Cercyonis Scudder, 1875; Pyronia 
Hiibner, [1819] and Aphantopus Wallengren, 1853 of Maniola Schrank, 1801 and Pseudonymphidia Callaghan, 1985 of 
Pachythone. Lafron Grishin, gen. n. (type species Papilio orus Stoll, [1780], parent subfamily Lycaeninae [Leach], [1815]) is 
described. Dipsas japonica Murray, 1875 is fixed as the type species of Neozephyrus Sibatani & Ito, 1942. The following taxa 
are junior subjective synonyms: Falcapica Klots, 1930 of Tetracharis Grote, 1898; Hahrodais Scudder, 1876, Favonius 
Sibatani & Ito, 1942, Neozephyrus Sibatani & Ito, 1942, Quercusia Verity, 1943, Chrysozephyrus Shirozu & Yamamoto, 1956, 
and Sibataniozephyrus Inomata, 1986 of Hypaurotis Scudder, 1876; Plesioarida Trujano & Garcia, 2018 of Roeberella Strand, 
1932; Papilio temenes Godart, 1819 (lectotype designated herein) of Heraclides aristodemus (Esper, 1794), Speyeria hydaspe 
conquista dos Passos & Grey, 1945 of Argymnis hesperis tetonia (dos Passos & Grey, 1945), and Erycides imbreus Plotz, 1879 
of Phocides polybius polybius (Fabricius, 1793). The following are revised genus-species combinations: Pachythone lencates 
(Hewitson, 1875) Pachythone flocculus (Brevignon & Gallard, 1993), Pachythone floccus (Brevignon, 2013), Pachythone 
heberti (P. Jauffret & J. Jauffret, 2007), Pachythone marajoara (P. Jauffret & J. Jauffret, 2007) and Cissia cleophes (Godman 
& Salvin, 1889). The following species are transferred between subgenera: Anthocharis lanceolata Lucas, 1852 belongs to 
Anthocharis Boisduval, Rambur, [Dumenil] & Graslin, [1833] instead of Paramidea Kuznetsov, 1929 and Danaus eresimus 
(Cramer, 1777) belongs to Danaus Kluk, 1780, and not to Anosia Hiibner, 1816. The following taxa are distinct species rather 
than subspecies (of species shown in parenthesis): Heraclides ponceana (Schaus, 1911) (not Heraclides aristodemus (Esper, 
1794)), Colias elis Strecker, 1885 (not Colias meadii W. H. Edwards, 1871), Argymnis irene Boisduval, 1869 and Argymnis 
nausicaa W. H. Edwards, 1874 (not Argynnis hesperis W. H. Edwards, 1864), Coenonympha California Westwood, [1851] (not 
Coenonympha tullia (Muller, 1764)), Dione incarnata N. Riley, 1926 (not Dione vanillae (Linnaeus, 1758)), Chlosyne 
coronado (M. Smith & Brock, 1988) (not Chlosyne fulvia (W. H. Edwards, 1879)), Chlosyne chinatiensis (Tinkham, 1944) 
(not Chlosyne theona (Menetries, 1855)), Phocides lilea (Reakirt, [1867]) (not Phocides polybius (Fabricius, 1793)), 
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Cecropterus nevada (Scudder, 1872) and Cecropteriis dobra (Evans, 1952) (not Cecropteriis mexicana (Herrich-Schaffer, 
1869)), Telegonus anausis Godman & Salvin, 1896, (not Telegonus anaphus (Cramer, 1777)), Epargyreus huachuca Dixon, 
1955 (not Epargyreus clarus (Cramer, 1775)), Nisoniades bromias (Godman & Salvin, 1894) (not Nisoniades rubescens 
(Moschler, 1877)), Pholisora crestar J. Scott & Davenport, 2017 (not Pholisora catullus (Fabricius, 1793)), Carterocephalus 
mandan (W. H. Edwards, 1863) and Carterocephalus skada (W. H. Edwards, 1870) (not Carterocephalus palaemon (Pallas, 
1771)), Amblyscirtes arizonae H. Freeman, 1993 (not Amblyscirtes elissa Godman, 1900), and Megathymus violae D. Stallings 
& Turner, 1956 (not Megathymus ursus Poling, 1902). Resulting from these changes, the following are revised species- 
subspecies combinations: Heraclides ponceana bjorndalae (Clench, 1979), Heraclides ponceana majasi L. Miller, 1987, 
Argynnis irene dodgei Gunder, 1931, Argynnis irene cottlei J. A. Comstock, 1925, Argynnis irene hanseni (J. Emmel, T. 
Emmel & Mattoon, 1998), Argynnis nausicaa elko (Austin, 1984), Argynnis nausicaa greyi (Moeck, 1950), Argynnis nausicaa 
viola (dos Passos & Grey, 1945), Argynnis nausicaa tetonia (dos Passos & Grey, 1945), Argynnis nausicaa chitone W. H. 
Edwards, 1879, Argynnis nausicaa schellbachi (Garth, 1949), Argynnis nausicaa electa W. H. Edwards, 1878, Argynnis 
nausicaa dorothea (Moeck, 1947), and Argynnis nausicaa capitanensis (R. Holland, 1988), Argymnis zerene atossa W. H. 
Edwards, 1890, Dione incarnata nigrior (Michener, 1942), Chlosyne coronado pariaensis (M. Smith & Brock, 1988), 
Cecropterus nevada aemilea (Skinner, 1893), Cecropterus nevada blanca (J. Scott, 1981), Telegonus anausis annetta (Evans, 
1952), Telegonus anausis anoma (Evans, 1952), Telegonus anausis aniza (Evans, 1952), Epargyreus huachuca profugus 
Austin, 1998, Carterocephalus mandan mesapano (Scudder, 1868) and Carterocephalus skada magnus Mattoon & Tilden, 
1998. American Coenonympha subspecies placed under C. tullia other than Coenonympha tullia kodiakW. H. Edwards, 1869, 
Coenonympha tullia mixturata Alpheraky, 1897 and Coenonympha tulliayukonensis W. Holland, 1900 belong to C. California. 
Heraclides ponceana latefasciatus Grishin, ssp. n. is described from Cuba. Argynnis coronis carolae dos Passos & Grey, 1942 
is considered a subspecies-level taxon. Unless stated otherwise, all subgenera, species, subspecies and synonyms of mentioned 
genera and species are transferred together with their parent taxa, and others remain as previously classified. 

Key words: taxonomy, classification, genomics, phylogeny, biodiversity. 
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INTRODUCTION 

DNA-based phylogenies (Zuckerkandl and Pauling 1965) revolutionized the way we view animal 
classification and taxonomy, including butterflies (Mutanen et al. 2010). Previously relying on several 
carefully selected gene markers, DNA methods have evolved towards sequencing and comparison of 
whole genomes. Genome-scale approaches aim at utilizing all DNA of an organism and thus are most 
comprehensive and accurate, frequently revealing inconsistencies between phylogeny and current 
classification (Kawahara and Breinholt 2014; Espeland et al. 2018; Allio et al. 2019; Li et al. 2019; Zhang 
et al. 2019a; Zhang et al. 2019b). Our research group genome-sequenced all butterfly species recorded 
from the United States and Canada (USC) (Zhang et al. 2019d) and proposed refinements to butterfly 
taxonomy (Zhang et al. 2019c). Currently we are working on extending our genomic datasets to cover 
subspecies and populations, in addition to species from other parts of the world. Review of these datasets 
suggests additional changes based on the analysis of phylogenetic trees and genetic variation between and 
within species. Here, we present these results after explaining the general logic and methods behind them. 

Total DNA is extracted from a leg (or other parts of a butterfly), fragmented into short pieces 
(unless the specimen is old and DNA is already degraded) and sequenced in 150 base pair (bp) segments. 
Thus, every fragment of DNA present in the sample is getting sequenced, resulting in the coverage of the 
entire genome, both nuclear and mitochondrial, coding and non-coding (Zhang et al. 2019a). Because this 
approach can sequence very short DNA (e.g. 25-50 bp), even very old specimens with degraded DNA 
yield usable data (Cong et al. 2019a). In this work, we mostly use DNA sequences of protein-coding 
genes, which are computationally selected from the total pool of DNA sequences to match a known 
protein set of a butterfly. Butterflies typically have nearly 15 thousand genes, and the total number of base 
pairs that we get for phylogenetic analysis is about 10 million (Cong et al. 2015; Zhang et al. 2019d). For 
comparison, a good dataset of gene markers typically covers less than 10 thousand base pairs. Due to the 
vast size of genomic datasets, phylogenies resulting from them are usually reliable. Reliability of each 
node in the tree is indicated by a number next to that node, the closer it is to 1, the more reliable is the 
node. Furthermore, the genotype encodes phenotype, thus, all the features of wing patterns and 
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morphology of butterflies and their life stages (including caterpillar foodplant preferences) are encoded in 
these genomic sequences. For this reason, the genomic dataset we sequence may be the best 
representation of a butterfly (richer than a pinned adult specimen that we visually inspect), given a 
meaningful approach to analyze these genomes. 

Armed with genome-scale phylogenetic trees, we have ventured to propose taxonomic changes. 
The majority of the changes discussed below deal with a change of rank in the classification. For genus- 
group names, we have changed the rank of taxa from genus to subgenus, or move between valid names 
and synonyms. For species-group names, mainly we have changed the ra nk of subspecies to species. The 
reason for these changes is to bring these categories in better agreement with their definitions, although 
such definitions are never precise and are continuously being refined with additional knowledge gained. 

A genus is defined as a monophyletic group (a clade in an accurate phylogenetic tree) of related 
species that is below the level of a tribe or subtribe. Currently, there are no objective criteria for grouping 
species into genera. However, the genus level is probably more important than any other classification 
level above species, because a species name contains a genus name as the first word. Therefore, for its 
optimal practical use, a genus should be defined neither too stringently to encompass only a small group 
of very close relatives, nor should it be too broadly defined to include too distant relatives. Subjectively, 
taxonomists follow their own intuition about where to draw a boundary of a genus. Our thoughts on the 
criteria for a genus are discussed in the Taxonomic Appendix to Li et al. (2019) and the Introduction in 
Grishin (2019). First, a genus should be a major and prominent phylogenetic group (i.e., it is best if the 
phylogenetic tree branch leading to the last common ancestor of the genus is longer than the nearby 
branches). Secondly, genetic diversity is a function of the time since the origin of a genus. Theoretically, a 
genus may be defined as all ancestors of a species that lived at some given point in the past. Practically 
this definition can be applied through a narrow slice across a time-calibrated phylogenetic tree at a 
specific time point (Talavera et al. 2012; Li et al. 2019). While it is not possible to establish a distance 
cutoff using COI barcodes (or any other gene markers) when defining a genus due to the reasons 
explained by Trujano-Ortega et al. (2020), a comparison of COI distances may be instructive and should 
be used as evidence combined with other considerations. Provided these general criteria, it is also best if 
most genera are in agreement with how they are currently defined (mostly using phenotypic 
considerations, such as similarity in appearance), to avoid additional name changes. Further examples of 
how we apply these principles to specific cases are given in the Discussion section below. 

Meaningful groups of species within a genus can be given a rank of subgenus. Subgenera are 
useful to define clades that are not as prominent as a genus, but represent groups of very close relatives, 
mostly recognizable immediately by their appearance as such. Although frequently frowned upon and 
synonymized or treated as a genus instead, a subgenus is one of just eight levels of the ICZN Code 
hierarchy of names (ICZN 1999). The number of levels in a phylogenetic tree is larger than these eight 
categories. Therefore, it seems wasteful to ignore the level that is aimed at refining the classification and 
indicating phylogenetic substructure within a large and diverse genus. We use subgenus category freely 
and suggest that many groups currently viewed as genera may be subgenera instead: it seems more 
valuable to stress their relatedness by uniting them in a broader genus, yet indicating their distinction by 
keeping them as subgenera. We believe that the consistency criterion is of importance here as well. I.e., 
groups defined as subgenera in one lineage should correspond to subgenera of comparable genetic 
differentiation and time of origin in another lineage. 

Species definitions (=concepts) have been extensively addressed in the literature (Mallet 1995). 
Comparing many species concepts that have been proposed (Aldhebiani 2018), our view is closest to what 
Claridge (2017) described as a "broadly based biological species concept", which seems similar to the 
"genomic integrity species definition" of Sperling (2003). I.e., species are defined by some reproductive 
barrier, which is not absolute but porous (Mallet et al. 2016); hybrids are characterised by lower fitness 
and usually are eliminated from the population, thus allowing each species to maintain its genetic 
uniqueness (integrity) that persists in time despite on-going hybridization with other species. This 
definition does not include phenotypic distinction between species, and allows for a possibility of 
"cryptic" species that look superficially indistinguishable, but are to some extent reproductively isolated 
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from each other. Equally, this definition does not exclude hybrid species, i.e. those that originate by a 
significant gene influx of one species into a population of another that occurred over a relatively short 
period of time. If this population that experienced the influx persists in time and is reproductively isolated 
from both of its parent species to the extent comparable to isolation in other related species, then it can be 
called a hybrid species, in particular, if it expands its range and further diversifies. 

Traditionally, species are defined by morphological distinctiveness. Morphological distinctiveness 
accumulates with time as a result of reproductive isolation, but it does not necessarily contribute to it. 
Therefore, using a morphological consideration does not provide direct evidence of reproductive isolation 
and thus speciation. On the contrary, using genomic DNA sequences, we can study reproductive isolation 
directly and look for segments of DNA from different species in a species’ genome under study. Presence 
of such segments would indicate hybridization and the fraction of such segments would indicate the 
extent of hybridization. Thus, analysis of genomic sequences may be the best practical approach to 
species definition available today. 

We carried out a study to devise genomic criteria for distinguishing different species from 
populations of the same species (Cong et al. 2019b). We studied sister butterfly populations across a 
central Texas suture zone. Suture zones are geographic boundaries common for many pairs of closely 
related species (counterparts) (Remington 1968). This central Texas suture zone is old (Newton 2003), 
therefore pairs that attained species status across it are likely to be more strongly isolated reproductively 
than more recently speciated pairs. Therefore, our criteria obtained using this old zone are likely to be 
conservative. We analyzed transcriptomes (mostly protein-coding RNA) of 25 pairs of western and 
eastern counterparts around the central Texas suture zone. Some of these counterparts were different 
species lik e Heraclides cresphontes (Cramer, 1777) and H. rumiko Shiraiwa & Grishin, 2014, others were 
populations of the same species from the north and south Texas, e.g. Hylephila phyleus (Drury, 1773). 

We found that two DNA-based measures best separate distinct species from conspecific 
populations. These measures were fixation index (Fst) and the extent of gene flow (measured by Gmin), 
both computed on the sex Z chromosome-encoded genes. Male butterflies have two Z chromosomes, 
while females have one Z and one W, which is a special female chromosome. Fst is frequently used in 
population genomics. Fst compares DNA divergence of specimens within each population (or species) to 
DNA divergence between populations (species). Fst ranges between 0 and 1. If Fst between two 
populations is low, near 0, then these populations are similar to each other. If Fst is high, above 0.5, then 
these two populations differ and are likely to be different species. Gmin is the fraction of segments of a 
different population (species) in a genome of a given population (species), thus it gives an estimate of 
gene exchange between populations or species (Geneva et al. 2015). Values of Gmin near 0 indicate the 
lack of gene exchange between populations suggesting that these populations are different species. Higher 
Gmin corresponds to more gene exchange and the lack of reproductive barrier meaning conspecificity. 

These two measures (Fst and Gmin) partition the 25 pairs of counterparts into two distinct and 
well-separated groups: the one corresponding to little genetic isolation and thus being the same species, 
and the other characterized by higher Fst (more than 0.2) and smaller Gmin (less than 0.06, i.e. about 6%), 
which corresponds to distinct species. Studies on human populations show that maximal Fst between 
human populations is about 0.2 (Nelis et al. 2009) and our species {Homo sapiens) has about 1.5% - 2.1% 
genes introgressed from Neanderthals {Homo neanderthalensis ) (Wall and Yoshihara Caldeira Brandt 
2016). Therefore, our criteria applied to humans confirm that modem humans are the same species, but 
Neanderthal was a distinct species. 

In this work, we compute Fst/Gmin on population pairs that we test for a possibility that they 
represent distinct species. When the values are, for instance 0.45/0.01, then it is most likely that we are 
dealing with distinct species (Fst is above 0.2 and Gmin is below 0.06). However, if the values we get are 
0.16/0.2, then these populations are conspecific (Fst is below 0.2 and Gmin is above 0.06). These criteria 
are useful, but they should not be taken absolutely and separately from other evidence. While large Fst 
values above 0.5 and small Gmin values below 0.01 are more definitive indicators of speciation, when 
they are closer to the "gray zone" (e.g. Fst is between 0.18 and 0.25 and Gmin between 0.03 and 0.07) 
further evidence is necessary. Additional problem may arise with siblings that are low in genetic diversity, 
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because their genomes are various recombinations of their parents' genomes, or with inbred populations, 
for a similar reason. Inclusion of these closely related individuals elevates Fst due to low diversity within 
siblings or inbred populations. To avoid this undesirable effect, we used specimens from different 
localities when possible. 

We also analyzed COI barcode divergence, which for distinct species is typically above 2% in the 
presence of phenotypic distinctiveness (Hebert et al. 2003) and inspect genomic trees. It should be taken 
into account that less than 2% barcode divergence has been reported for distinct species (Burns et al. 
2007). Conversely, due to introgression, conspecific individuals may differ in their barcodes by more than 
2% (Zakharov et al. 2009; Cong et al. 2017). 

Furthermore, distinct species usually form distinct and well-supported clades (support values near 
1 by these clades) in phylogenetic trees. Branches supporting distinct species are typically longer than 
internal branches within a subtree of conspecific specimens. Each terminal branch in the tree is usually 
long, because of individual variation and uniqueness of each specimen, and also because sequencing 
errors and contaminations (not being common for any pair of specimens) tend to be reflected in the 
terminal branch. For conspecific populations, specimens may be intermixed in the trees, or support values 
for separation of populations are smaller (closer to 0 than to 1), due to significant gene exchange between 
these populations, when some genes cluster individuals differently from other genes. 

Yet another advantage of genomic analysis is its robustness to small number of specimens 
analyzed. Phenotypic analysis must rely on a large number of specimens to gauge the range of variation, 
which is shaped by the interactions of genotype with the environment. DNA variation is not affected by 
the enviromnent, and each specimen contains two genomes: from mother and from father, which in turn 
contain information from "grandparents". Therefore, even a couple of specimens from distant localities 
(i.e. no close kinship) is sufficient to estimate intraspecific variation with reasonable accuracy. 

For brevity, we do not review phenotypic differences between most of the taxa we deal with, and 
interested readers should consult other publications. All these taxa have been previously described and 
their phenotypic characters given elsewhere. These characters do not contribute to our decision to change 
ranks of these taxa, although they may be looked at as complementary supporting evidence that increases 
confidence in the results. The evidence presented here is based on genomic analysis. The data we offer are 
new and they allow us to interpret known phenotypic differences between these taxa from a different and 
complementary perspective. The purpose of this work is to propose taxonomic changes that are gleaned 
from genomic data, so that the new name combinations can be used in other publications. A more detailed 
evolutionary analysis of genomic data will be published elsewhere. 

The taxonomic rearrangements we propose are supported by genome-scale phylogenetic trees, Fst 
and Gmin statistics, and sometimes complemented by phenotypic considerations. The sections are 
presented in standardized format. The taxonomic act is the title of each section. Relevant genera, 
subgenera and their type species, are specified. When the species are listed with their original genus 
name, author names are given without parenthesis. For each species and subspecies with revised rank, 
type locality is given. A section is usually illustrated with a small segment of a nuclear genomic tree (or 
other trees as stated in the text) including taxa necessary to support the conclusions. Previous, not newly 
proposed, genus-species and species-subspecies combinations are used in the figures (per Pelham 2020 
<http://www.butterfliesofamerica.com/US-Can-Cat.htm>, version revised 7 August 2020, except those 
name changes adopted before this publication based on our genomic results). New name combinations are 
given in the text. Taxa of major focus are shown in red, other taxa of interest are shown in blue, magenta 
or green. The section ends with a conclusion and in many cases with a list of species with revised genus- 
species or species-subspecies name combinations. The sections are ordered by family and typically in 
their taxonomic order from Zhang et al. (2019d) (except in Nymphalidae, where the arrangement was 
altered to better fit larger images on pages while keeping them next to relevant text). Finally, whole 
genome shotgun datasets we obtained and used in this work are available from the NCBI database 
<https://www.ncbi.nhn.nih.gov/> as BioProject PRJNA672791, and BioSample entries of the project 
contain the locality and collection data of the specimens sequenced. 
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Family Papilionidae Latreille, [1802] 

Heraclides ponceana (Schaus, 1911) is a species 
distinct from Heraclides aristodemus (Esper, 1794) 

Initially proposed as a species and currently treated as a subspecies of Heraclides aristodemus (Esper, 
1794) (type locality Haiti), Papilio ponceana Schaus, 1911 (type locality USA: Florida, Miami) shows 
profound genetic differentiation from H. aristodemus by a magnitude characteristic of species-level taxa. 
Phylogenetic trees reveal partitioning of the taxa previously placed in H. aristodemus into two prominent 
clades, including the tree constructed from proteins encoded by the Z chromosome (Fig. 1). The Fst 
between H. a. ponceana and H. 
aristodemus populations is 0.43, 
which indicates their strong 
differentiation, and Gmin is 0.02, 
which suggests low gene exchange 
between them. COI barcodes of the 
two also differ by 3.2% (~21 
different base pairs). Moreover, it is possible (see below) that H. a. ponceana and H. aristodemus are (or 
were) sympatric in Cuba. Therefore, we reinstate Heraclides ponceana (Schaus, 1911) as a species-level 
taxon. Due to their genetic and morphological similarities, we consider Heraclides ponceana bjorndalae 
(Clench, 1979) (type locality Bahamas: Great Inagua Island) and Heraclides ponceana majasi F. Miller, 
1987 (type locality Bahamas: Crooked Island) to be subspecies of H. ponceana, new combinations. 

Papilio lemenes Godart, 1819 is a junior subjective synonym 
of Heraclides aristodemus (Esper, 1794) 

The short version of the original description of Papilio temenes Godart, 1819 (page 21) is: "Dessus des 
ailes d’un brun-noiratre, avec deux bandes jaunes, maculaires, disposees en sautoir sur les superieures: les 
inferieures en queue: le dessous de celles-ci jaunatre, avec une bande bleue, flexueuse, sur le milieu" 
(Godart 1819). It can be translated word-for-word as: "Above the wings of a blackish-brown, with two 
bands yellow, macular, arranged in sautoir [=crosswise] on the forewings: the hindwings [end] in tail: the 
underside of these [hindwings, not tails] yellowish, with a band blue, flexuous, in the middle", and 
interpreted as: "Wings above blackish-brown, forewings with two yellow macular bands that cross each 
other, hindwings with tails, underside of hindwings yellowish with a blue flexuous band in the middle." 
Perhaps the most significant word here is "sautoir" defined in Dictionary.com as "A ribbon, chain, scarf, 
or the like, tied around the neck in such a manner that the ends cross over each other" (Dictionary.com 
FFC 2020), and by Merriam-Webster dictionary as "a chain, ribbon, or scarf worn about the neck with the 
ends forming a St. Andrew's cross in front" (Merriam-Webster 2020). It is important because this 
character (crosswise vs. more parallel arrangement of yellow forewing bands) is diagnostic of H. 
aristodemus (crosswise) vs. H. ponceana (more parallel, images in Warren et al. (2016)). Today, the name 
temenes is applied to a broad-banded form of H. ponceana that flies in Cuba (see below). Here we argue 
that this is a misidentification, and the Godart’s temenes is H. aristodemus instead. 

We see that, although the Godart's description of P. temenes is brief, two of its details agree better 
with H. aristodemus than with the broad-banded subspecies of H. ponceana from Cuba: crisscrossing 
bands and blue flexuous band. In the broad-banded subspecies, the outer band is more parallel to the inner 
band and does not give an impression of crossing it. And the blue spots on the hindwing underside look 
more like a row of lunules than an irregular blue band. In H. aristodemus, forewing bands indeed give an 
impression of crossing bands, and the blue spots on the hindwing below look like an irregular blue band. 
Furthermore, an extended description of P. temenes on page 63 states that the forewing bands are narrow 
and macular ("etroites, maculaires") (Godart 1819), instead of being broad and continuous. Besides this 
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additional detail, the extended description reiterates other points of the short description. We note that 
some H. aristodemus females may have broader and more continuous yellow bands that are more similar 
to those in H. ponceana (Fig. 2e), but the P. temenes description does not mention this possibility, 
simplifying the application of the name. 

Next, we inspected all potential syntypes of P. temenes, two specimens in Paris, France (Fig. 2ab) 
and one (only photograph inspected) in Edinburg, UK (Bland 2019). These 3 specimens closely agree 
with the original description of P. temenes, carry historical labels and labels indicating their type status 
("type", "co-type", or "?co-type") and therefore are likely to be true syntypes of this taxon. These 
specimens phenotypically are H. aristodemus and not H. ponceana, differing from the broad-banded 
subspecies from Cuba. Furthermore, we sequenced two syntypes in Paris (one labeled "TYPE", the other 
labeled "CO-TYPE") and they are H. aristodemus by genomic DNA (Fig. 1 blue), in agreement with their 
wing patterns. In conclusion, our analysis of the original description and the likely syntypes leaves little 
doubt about the identity of P. temenes as II. aristodemus, which is a species different from the broad- 
banded H. ponceana "temenes " found on Cuba today (Fig. 1 red). 




Fig. 2. Heraclides aristodemus (a-e) and II. ponceana (f-h). Green arrows point at equivalent positions in different specimens. 


Due to misapplication of the name temenes to a taxon different from that in the original description 
and for the future stability of the name, we designate the syntype in the Museum National d'Histoire 
Naturelle, Paris, France (MNHP) shown in Fig. 2a and possessing all the characters stated in the original 
description, with the following 4 rectangular white (some faded to brownish) labels: printed red || TYPE ||; 
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framed, lined and printed with green, first line handwritten in black || Anc. Collection | MUSEUM DE 
PARIS ||; printed, with a square barcode on the right side || MNHN, Paris | EL63126 ||; printed || DNA 
sample ID: | NVG-18079F07 | c/o Nick V. Grishin ||, the lectotype of Papilio temenes Godart, 1819. The 
red, rectangular, printed label || LECTOTYPE I Papilio temenes | Godart, 1819 | designated by Grishin 
|| will be added to this specimen. This specimen was chosen as the lectotype because it is labeled "type" 
rather than "co-type", and it yielded genomic dataset of a good quality for a specimen that old. 

It appears that Oberthiir (1897) was the first to incorrectly apply the name P. temenes to the broad- 
banded subspecies of 77. ponceana from Cuba, inconsistently with the original description and now with 
the identity of the lectotype of P. temenes. Oberthiir illustrated one such specimen. While we have not 
investigated the reasons behind Oberthur's mistake, we note that it has been followed in subsequent 
literature. Interestingly, we found 3 century-old 77. aristodemus (i.e. true P. temenes) specimens labeled 
"Cuba", in ZMHB and FMNH (Fig. 2c-e). Further studies may answer the question whether both species 
(77. aristodemus and 77. ponceana) co-occurred in Cuba, or the old records from Cuba were mislabeled. 

In summary, we conclude that Papilio temenes Godart, 1819, syn. n., is a junior subjective synonym 
of Heraclides aristodemus (Esper, 1794). The type locality of temenes should remain as stated on page 63 
(Godart 1819): "Antilles & dans l’Amerique septentrionale", i.e., "West Indies and North America", which 
is not necessarily Cuba. It remains to be investigated whether 77. aristodemus has been found or still 
occurs in Cuba. In any case, as detailed above, the broad-banded subspecies of 77. ponceana from Cuba 
does not have a name, and it is described here as new. 

Heraclides ponceana latefasciatus Grishin, new subspecies 

http://zoobank.org/C835D721-548C-4822-9280-E241AlC94866 

(Figs. 1, 2f-h) 

Definition. This taxon differs from 77. aristodemus in that the outer band of yellow spots on the forewing 
is more parallel to the inner band (and the outer margin) rather than approaching the inner band at an 
angle and giving an appearance of crossing the inner band; on the hindwing below there are more 
prominent red spots and more crescent-shaped blue spots that are better separated from each other, rather 
than fonning a continuous band. This subspecies differs from all other 77. ponceana subspecies by broader 
yellow central bands on both wings, less extensive brown coloration on the forewings below, a paler basal 
area and the lack of red spotting in the postdiscal area on hindwing above. 

Type locality. Cuba: Guantanamo province, Rio Seco, San Carlos Estate. 

Distribution. Known only from Cuba. 

Etymology. The broad yellow bands are the most distinctive feature of this subspecies. The name is 
formed from the Latin words latus (wide, broad) and fascia (band, stripe). The name is an adjective. 

Type material. Holotype male (Fig. 2g), deposited in the American Museum of Natural History, New 
York, NY, USA (AMNH), with the following 4 rectangular white (some faded to brownish) labels: 
printed, the date crossed out || San Carlos Est. | Guantanamo | Cuba. 4 8 X T3 ||; handwritten || May 1 '00 
||; printed with the numbers handwritten || Am. Mus. Nat. Hist. | Dept. Invert. Zool. | No. 20976 ||; printed 
|| DNA sample ID: | NVG-14101H09 | c/o Nick V. Grishin ||. The red, rectangular, printed label || 
HOLOTYPE S I Heraclides ponceana \ latefasciatus Grishin || will be added to this specimen. Ten 
paratypes (when known, localities are given in parenthesis after specimen numbers): 5 males (NVG- 
14101H10 in AMNH and NVG-14106A11 (Matanzas), NVG-14106A12 (Santiago), NVG-14106B01 
(Guantanamo) & NVG-14106B02 in USNM) and 5 females (NVG-14114B09 & NVG-14114B10 in 
LACM and NVG-14106B03 & NVG-15104A01 (both from Santiago) & NVG-15104A02 in USNM). 

Barcode sequence of the holotype. aacattatattttatttttggtgtttgagcaagaatattaggaacttctcttagtttatta 

ATTCGAACTGAATTAGGAACTCCAGGTTCTTTAATTGGAGATGATCAAATTTATAATACCATTGTTACAGCTCATGCTTTTATTATAATTTTTT 
TTATGGTTATACCTATTATAATTGGAGGATTTGGTAATTGATTAGTTCCATTAATATTAGGAGCCCCTGATATAGCTTTCCCTCGAATAAATAA 
TATAAGATTTTGACTTTTACCTCCTTCTTTAACTCTTTTAATTTCAAGTATAATTGTCGAAAATGGAGCTGGAACTGGATGAACTGTTTATCCT 
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CCCCTTTCTTCTAATATTGCTCATGGAAGAAGTTCAGTAGATTTAGTTATTTTTTCTCTTCATTTAGCGGGTATTTCTTCAATTTTAGGAGCAA 
TTAATTTTATTACTACTATTATTAACATGCGAATTAATAGAATATCCTTTGATCAAATACCTTTATTTGTTTGAGCTGTAGGAATTACAGCTTT 
ATTATTACTCTTATCCTTACCCGTTTTAGCTGGAGCTATTACTATATTATTAACTGATCGAAATTTAAATACTTCATTCTTTGATCCTGCAGGA 
GGAGGAGATCCTATTCTATACCAACACTTATTT 

Instead of proposing a new name for the Cuban broad-banded subspecies of H. ponceana, we entertained 
a possibility to request ICZN to designate one such specimen as the neotype of P. temenes, consistent 
with the current usage of this name, but contrary to the original description and the identity of three extant 
syntypes. However, we decided against this route for the following reasons. First, the name H. a. temenes 
is not in very wide use being applied to an uncommon endemic of a single island to warrant a special 
consideration by ICZN. Second, it seems most fair to respect original research that lead to creation of this 
name, and the original identity of this species that is quite clear even from its description alone (see 
above). Third, it is conceivable that the true P. temenes occurred (or even still occurs) in Cuba, and further 
research may show that it is not a synonym, but a valid subspecies of H. aristodemus, in which case it will 
be without a name if the neotype is designated to preserve the current usage of "temenes" as a subspecies 
of H. ponceana. It would create a nuisance situation when the original P. temenes would need a new 
name. Finally, a valid name that is suggestive of a diagnostic character ( latefasciatus ) has an advantage of 
being easier to attribute to the taxon (compared to temenes ) and thus may be easier to remember. 


Family Pieridae Swainson, 1820 

Colias elis Strecker, 1885 is a species 
distinct from Colias meadii W. H. Edwards, 1871 

Previously considered a subspecies of Colias meadii W. H. Edwards, 1871 (type locality USA: Colorado), 
Colias elis Strecker, 1885 (type locality 
Canada: Alberta) is not monophyletic with 
C. meadii in the genomic trees (Fig. 3, 
although with low support), including the 
tree constructed from Z chromosome- 
encoded genes. The identity of the C. elis 
clade is supported by the sequence of its 
syntype from the Field Museum National 
History collection (Fig. 3). Furthermore, C. 
elis and C. meadii show genetic differences 
that are larger than a number of other 
Colias species pairs, e.g. Colias eriphyle W. 

H. Edwards, 1876 (type locality Canada: 

British Columbia) and Colias eurvtheme 
Boisduval, 1852 (type locality USA: 

California). Fst/Gmin statistics for the comparison of C. elis and C. meadii are 0.40/0.03, but those for C. 
meadii meadii and C. meadii lemhiensis are 0.18/0.14 (can be used as a control of conspecific taxa). 
Therefore, we reinstate Colias elis Strecker, 1885 as a species-level taxon. 




Colias nastes streckeri|17114D02|M|Canada:Alberta|2003 
— Colias boothii boothii 17114C12 M|Canada:Nunavutl 1988 
- Colias heda|17114C09lM|USA:AK,Yukon-Koyukuk 1979 
• Colias rankinensisTl9035F10lM|Canada:NWT,Schultz Lake Rapids|1965 

L 244 —- - Colias johansenil 19062H07lPT|F|Canada:NWT| 1988|LACM 

LTItt^ - Colias johanseni|17114C07|F|Canada:Northwest Territories! 1988 

0.52 - Colias hyperborea|19043G03|F|Russia, Siberia 

- Colias meadii elisll8039F07|ST|F|Canada:Alberta|prior to 1885|FMNH 

— Colias meadii elisJ19035Fll|M|Canada:Alberta,East Cirque|1977 
Colias meadii elis I l9047B08]MjCanada:AB,Kananaskis For. Res.|2001 
Colias meadii elis|19035F12|F|Canada:Alberta,Rocky Mtns., Jasper Pk.| 1986 
Colias meadii elis119047B09|F|Canada:AB,Prospect Mountain! 1989 
Colias meadii elis119063Fll|M|Canada:Britisn Columbial 1984 
as meadii lemhiensis|19063G05|MlUSA:WY,Park Coj2007 
Colias meadii lemhiensisl 19036E12]F|USA:ID,Lemhi Co., Lemhi Pass|1978 
Colias meadii lemhiensis|9394|M|USA:WY,Park Co^Yellowstone NP|2017 
Colias meadii meadii[PAOl59|USA:CO,Clear Creek Co., Loveland Pass|2016 
Colias meadii meadiill9063G04|M[USA:WY < Albany Co.|2018 
— Colias meadii meaaiiJ19063G06lM|USA:CO,Routt Co.. Muddy Pass11998 
Colias gigantea mayi|17114D04|M|Canada:Alberta|2003 
Colias scudderiilPAO160|USA:CO,Clear Creek Co.,Loveland Pass|2016 
Colias occidental^ chrysomelas|GK6j_USA:CA,Colusa Co.,Mill Valley Camp|2017 
Colias Christina krauthii|19039A05.‘. 



5IM|USA:SD,Black Hills| 1986 


Colias alexandraapachej 11560J m]US_A:AZ,A pache_Co.[2018 


Colias harfordii|17114C06]MlUSA:CA,San Diego Co.J2008 
Colias vitabunda|18042E08|F[Canada:Yukon Territory12016 
Colias eriphyle194211MJUSA:V7Y.Park Co.,Yellowstone National Park|2017 
Colias philoaice14261 [MlUSA:IN,Montgomery Co.|2015 
Colias eurytheme|6530|F|USA:CO.Park Co.,4.5 air mi SSE Florissant|2016 
Colias palaeno chippewall9034A10lM|Canada:YT| 1982 
Colias pelidne|171l9E03lM|Canada:Newfoundland|1975 
Colias interior dornfeldi 18062B11|M|USA:MT,Missoula Co.|1994 
Colias skinneri117114D06IF]USA:WY,Fremont Co.|2007 
Colias behriil 18038A05|M|USA:CA,Mono Co.|2009 

Fig. 3. Colias elis (red), C meadii (blue) and others w/ patch (magenta). 


Tetracharis Grote, 1898 is a valid subgenus that includes Anthocharis cethura C. 
Felder & R. Felder, 1865 (Muller, 1764), Anthocharis midea (Hiibner, [1809]), and 

Anthocharis limonea (A. Butler, 1871) 

The genomic tree of representative species of Anthocharis Boisduval, Rambur, [Dumenil] & Graslin, 
[1833] (type species Papilio cardamines Linnaeus, 1758) including all type species of available genus- 
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group names considered subgenera or junior subjective synonyms of Anthocharis revealed an unexpected 
but highly confident clade consisting of Anthocharis cethura C. & R. Felder, 1865 (type locality USA: 
California, Los Angeles Co.), Anthocharis 
midea (Hubner, [1809]) (type locality USA: 

Georgia, Wilmington Island) and Anthocharis 
limonea (A. Butler, 1871) (type locality 
Mexico) (Fig. 4 red). The clade is 
unexpected, because currently A. cethura is 
placed in the subgenus Anthocharis, while the 
other two species are placed in the subgenus 
Paramidea Kuznetsov, 1929 (type species: Anthocharis scolymus Butler, 1866). Curiously, neither of 
these species in the red clade are monophyletic with the type species of subgenera they are currently 
assigned to: A. (Anthocharis) cardamines is in the blue clade (Fig. 4) and A. (Paramidea) scolymus is in 
the green clade. Therefore, the current classification is incorrect. Out of genus-group names that are 
available for the red clade, Tetracharis Grote, 1898 (type species Anthocharis cethura C. & R. Felder, 
1865) is older than Falcapica Klots, 1930 (type species Papilio genutia Fabricius, 1793, which is a junior 
homonym of Papilio genutia Cramer, 1779, the oldest available name for this species is Mancipium 
midea Hubner, [1809]). As a result, we treat Tetracharis as a valid subgenus, new status, that includes 
three species: A. cethura, A. midea, and A. limonea, making Falcapica its junior subjective synonym. 



Anthocharis sara sara|5638|USA:CA,San Diego Co. 
Anthocharis iulia julia|1952|USA:CO,Boulder Co. 
Anthocharis tnoosa thoosa|1970|USA:AZ.Mohave Co. 
Anthocharis lanceolata lanceolata|PA063[USA:CA,Sierra Co. 
Anthocharis belia|19099H10|Morocco 
Anthocharis cardamines|: 


Anthocharis cethura 
Anthocharis midea i 


s|17115H10|M|Spain:Catalonia 


cethura16667[USA:CA.San Diego Co 
annicLae[6134|USA:WV,Pendlefcon C 
- — "-il f)\c ' -- 


Co. 


Anthocharis limonea 129491 Mexico: Mexico 
j —r Anthocharis scolymus! 19067G08|Japan 
1 — Anthocharis scolymus|19067G09|South Korea 

Euchloe tagis|19099H06|Spain 

- Euchloe ausonides coloradensis|6465|USA:CO,Grand Co. 

- Euchloe lotta|PA03911USA:NV.Elko Co. 

Fig. 4. Subgenera Anthocharis (blue), Tetracharis (red) 
and Paramidea (green) with Euchloe as outgroups. 


Anthocharis lanceolata Lucas, 1852 belongs to subgenus Anthocharis Boisduval, 
Rambur, [Dumenil] & Graslin, [1833] instead of Paramidea Kuznetsov, 1929 

Our genome-level phylogeny strongly supports the placement of Anthocharis lanceolata Lucas, 1852 
(type locality "Californie") as sister to the Anthocharis sara Lucas, 1852 (type locality "Califomie") 
species group, which belongs to the subgenus Anthocharis Boisduval, Rambur, [Dumenil] & Graslin, 
[1833] (type species Papilio cardamines Linnaeus, 1758) (Fig. 4, blue). The A. sara group taken together 
with A. lanceolata is sister to the subgenus Tetracharis (type species Anthocharis cethura C. & R. Felder, 
1865) and thus is not monophyletic with Anthocharis scolymus Butler, 1866, the type species of 
Paramidea Kuznetsov, 1929. In other words, Paramidea Kuznetsov, 1929 is sister to a clade formed by 
subgenera Anthocharis and Tetracharis. Therefore, A. lanceolata does not belong to Paramidea, but 
instead is in the subgenus Anthocharis. Apparently, the falcate shape of the forewing vs. rounded 
forewing apex is not a character that indicates phylogenetic groupings within the genus Anthocharis and 
has originated more than once. 


Family Lycaenidae [Leach], [1815] 

Tharsalea Scudder, 1876, Helleia Verity, 1943, Apangea Zhdanko, 1995 and 
Boldenaria Zhdanko, 1995 are genera distinct from Lycaena [Fabricius], 1807 

Discovering non-monophyly of Lycaena [Fabricius], 1807 (type species Papilio phlaeas Linnaeus, 1760) 
in a pioneering DNA-based phylogenetic analysis (van Dorp 2004) but stopping short of imminent 
taxonomic lumps or splits, de Jong and van Dorp (2006) concluded: "we propose that first the 
interrelationships as suggested by the present study are confirmed by further genetic markers." We 
accepted this challenge and used not some further markers, but all protein-coding genes in 34 Lycaeninae 
species, including the type species of 22 out of 28 available genus-group names (two of which share valid 
name of the type species with another genus-group name). The results confirm that Lycaena is not 
monophyletic (Fig. 5). Notably, Iophanus Draudt, 1920 (type and the only species Chrysophanus (?) 
pyrrhias Godman & Salvin, 1887) originates within Lycaena and is sister to most other American species 
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with strong support. This placement is unexpected due to the prominent phenotypic similarities between 
Iophanus and Melanolycaena Sibatani, 1974 (type species Melanolycaena altimontana Sibatani, 1974). 
Even using genitalic morphology, Iophanus 
pyrrhias (type locality Guatemala) was 
associated with the (largely) Palearctic clade 
(Fig. 5 blue) by Klots (1936) who wrote: 

"pyrrhias ... relation-ship is undoubtedly with 
the Palaearctic rather than with the Nearctic 
series, and is possibly rather ancient." As the 
genomic tree suggests, Klots was incorrect on 
both counts: I. pyrrhias is a relatively recent 
offshoot of the (largely) Nearctic clade (Fig. 5 
red). To the contrary, Melanolycaena is sister 
to Lycaena boldenarum White, 1862 (type 
locality New Zealand) and they form a clade 
that is sister to all other Fycaeninae we have 
sequenced, including Heliophorus Geyer, 

[1832] (type species Heliophorus belenus 
Geyer, [1832], considered to be a junior 
subjective synonym of Polvommatus epicles 
Godart, [1824]). Thus, we were impressed by 
the intuition of Sibatani (1974) who stated in the last sentence of his work: "the possibility is not 
completely ruled out that the New Guinean Melanolycaena and the coppers of New Zealand are 
monphyletic [sic!]." Indeed, Melanolycaena is sister to Boldenaria Zhdanko, 1995 (type species Lycaena 
boldenarum White, 1862) from New Zealand (Fig. 5), despite phenotypic dissimilarities. For these 
reasons, genomic phylogeny implies that the division of Fycaeninae into two sections (. Lycaena and 
Heliophorus) as suggested by Eliot (1973) was indeed tentative and needs to be revised because these 
sections are not monophyletic. 

Here, armed with genomic data for ~80% of the available genus-group names, we attempt such a 
revision. One solution to restore monophyly is to treat the whole subfamily as a single genus Lycaena that 
subsumes Melanolycaena and Heliophorus among others. This super-lumping solution may be in 
agreement with relatively low genetic diversification among all these species. Indeed, Fycaeninae 
experienced some of the slowest evolutionary rates among Fycaenidae as revealed by relatively shorter 
branches within the Fycaeninae clade compared to others (Zhang et al. 2019d). The difference in COI 
barcodes between distantly related Fycaeninae species ranges from about only 5% (35 bp, L. boldenarum 
and L. phlaeas, a difference common for closely related congeners) to about 8% (53 bp, Heliophorus sena 
Kollar, 1844 and Lycaena pang Oberthur, 1886, a difference typical for distantly related congeners). This 
high similarity in COI underscores the idea that it is meaningless to impose a strict cutoff on divergence 
values due to the differences in evolutionary rates in different lineages. However, even American species 
of Lycaena are estimated to have diverged over 20 million years ago (Zhang et al. 2019d), which is larger 
than for a typical diverse genus; i.e., it is larger than the divergence between Anthocharis and Euchloe 
Hubner, [1819] (Pieridae) and about the same as between Vanessa [Fabricius], 1807 and Nymphalis Kluk, 
1780 (Nymphalidae). Furthermore, Heliophorus has been traditionally maintained as a genus-level taxon. 
Therefore, we reject the super-lumping solution of a monotypic subfamily Fycaeninae. 

The opposite extreme would be to find a meaningful level closest to the leaves of the tree that defines 
genera. Ideally, there would be situations in the tree where many lineages diverge at about the same level 
(i.e. at the same distance from the root, meaning at about the same time in the past) and then stay as single 
lineages for some time (i.e. form longer branches). This rapid diversification immediately followed by a 
relative lack of further diversification creates a level in the classification, i.e. the tree looks more like a 
bush or a comb than a bifurcating tree at that point. Taking the (largely) Nearctic group (Fig. 5 red), we 
see exactly this situation: at its base, this group diversifies into five prominent clades, and then two of 
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Lycaena boldenaruml 19079Dll|New Zealand|1979 

Melanolycaena altimontanal 19079E02|Papua New Guinea|1973 

. . e fj c _.------ 


Heliophorus sena| 19079E03|Pakistan| 1976 



Heliophorus epicles|19079D 121Malaysia:Pahang11991 
eliophorus yunnanil 19079E05|China:Yunnan|1996 
Heliophorus gloria119079E01|Myanmar:Kachin| 1994 
Heliophorus tamu| 19079E04|Himalaya|old 

ycaena pang| 19093B12|China:Sichuan| 1902 

Lycaena helle119079D10|Russia:Siberia,Sayan Mountains! 1970 


I- Helioph 

1 1 1 — Heliopho 
l—f - ! Heliophc 
' - Hplinn 


j-j Lycaena phlaeas ursadental 17115H05|USA:WY,Park Co.11997 
Lycaena phlaeas phlaeas|PAOE01|France|2017 
q go Lycaena alciphron|19093Bll|Ukraine,Kiev|1969 

Lycaena cupreus lapidicola|PA0475|USA:CA,Mono Co.|2017 
Lycaena aoelidesl 17066G09|PT|Tajikistan| 19981 
Lycaena thersamonl 19079D07|Tajikistan:Pamir| 1988 
Lycaena tityrus|PAOE02|France]2017 
Lycaena virgaureae|PAOE09|Switzerland|2017 
Lycaena disparl 19079D06|Ukraine,Kiev| 1967 
Iophanus pyrrhias| 17066Gll|Guatemala| 1967 
Lycaena caspius| 19079E06|Iran,Schakuh|old 

- Lycaena hermes|6685|USA:CA,San Diego Co.|2011 
Lycaena arota virginiensis|PA0449|USA:CA,Plumas Co.|2017 

Lycaena rubidusjPA0126|USA:WY,Laramie Co.|2016 
n Lycaena xanthoides| 17114Dll|USA:CA ; Santa Clara Co.|2007 
_ Lycaena editha|PAO60|USA:CA,Sierra Co.|2016 

0 . 74 — Lycaena dione| 10599|USA:TX, Hemphill Co.|1977 

- Lycaena hyllus|PA0286|USA:C0,Adams Co.|2016 
Lycaena heteronea|9348|USA:WY,Parh Co.|2017 
Lycaena gorgon gorgon|PA0346|USA:CA,Santa Cruz Co.|2017 

o.9_4 - 1 Lycaena nivalis browni|6459|USA:CO,Grand Co.|2016 
_y— Lycaena epixanthe epixanthe|6686|Canada:Quebec|2011 
r 1 Lycaena mariposa penroseae|9440|USA:WY / Park Co.|2017 
Hr* Lycaena helloides|19077Bll|USA:CA,San Benito Co.11991 
IpoUfaena dospassosil 17114D12|Canada:New Brunswick|1983 
Lycaena dorcas dorcas|19081G05|Canada:Manitoba|1979 

Fig. 5. Lycaena (blue), Tharsalea (red), Iophanus (green), Helleia 
(magenta), Apangea, (cyan), Melanolycaena (purple), Heliophorus 
(orange-brown) and Boldenaria (black). 































these clades diversify further, also at approximately the same time point in the past. These five clades 
form a level in the tree and can be used as genera, offering the splitting solution. Notably, every one of 
these clades already has a genus-group name (Pelham 2008), including Palearctic Hyrcanana Bethune- 
Baker, 1914 (type species Polyommatus caspius Lederer, 1870). Apparently, these clades were also 
obvious from phenotypes: that is how they were defined and named to begin with (Scudder 1876; Klots 
1936; Miller and Brown 1979). This level of classification can be propagated to other parts of the tree, 
although they are currently poorly covered by species. It is a meaningful level that can be chosen to define 
genera, but a significant number of such genera will be monotypic (e.g. two out of five in the Nearctic 
clade), and excepting knowledgeable aficionados of this group, such genera carry little information about 
their interrelationships. Hence, we looked for a compromise between the splitting and lumping solutions. 

Inspection of the tree reveals a rapid diversification point between its root and the diversification of 
the red clade (Fig. 5): i.e. orange-brown, cyan, magenta, blue and red + green clades diverged at 
approximately the same time in the past. This divergence is followed by the lack of immediate further 
divergence, creating long and prominent branches in the tree and resulting in a meaningful level for 
classification. We have chosen to take this intermediate level to suggest division of Lycaeninae into 
genera. Most of these genera are unambiguously apparent from the tree: black, purple, orange-brown, 
cyan, magenta, blue and red clades stand for seven genera, all of which have previously proposed names. 
Two instances require further elaboration. First, Heliophorus sena (Kollar, [1844]) stands out from the 
rest in the genus (Fig. 5 orange). The type species of subgenus Nesa Zhdanko, 1995, it may be a genus- 
level taxon. However, it is monophyletic with Heliophorus and we leave it there as a subgenus to 
emphasize this relationship, awaiting further studies. Second, Iophanus (Fig. 5 green) is at about the same 
divergence from the rest of Nearctic species (Fig. 5 red) as H. sena from other Heliophorus. For now, we 
decided to keep this monotypic genus, because it is currently treated as such, and because its earlier 
divergence time sets it apart from the rapid diversification of the red clade. The name for the red clade is 
Tharsalea Scudder, 1876 (type species Polyommatus arota Boisduval, 1852), as chosen by Klots (1936), 
probably because this name was proposed before others in the paper (Scudder 1876). 

In summary, we refrain from partitioning Lycaeninae into tribes and revise the status of the 
following names treating them as genera: Tharsalea Scudder, 1876, Helleia Verity, 1943 (type species 
Papilio helle Denis & Schiffermiiller, 1775), Apangea Zhdanko, 1995 (type species Chysophanus pang 
Oberthur, 1886) and Boldenaria Zhdanko, 1995. Furthermore, in agreement with previous studies 
(Sibatani 1974; van Dorp 2004; de Jong and van Dorp 2006), we conclude that the endemic South African 
species currently placed in Lycaena represent the 9th genus of Lycaeninae that is named next. We are 
looking forward to testing this hypothesis with genomic data. 


Lafron Grishin, new genus 

http://zoobank.org/0DB9D8C5-E666-46A3-822B-50E96448C82A 
Type species. Papilio orus Stoll, [1780]. 

Definition. In male genitalia (Fig. 2 in de Jong and van Dorp 2006), differs from others in the subfamily 
Lycaeninae, except Melanolycaena, by a saccus-like pouch on juxta (Sibatani 1974); separated from 
Melanolycaena (Fig. 4 in Sibatani 1974) by juxta connected to valva at its more ventral part, as in other 
Lycaena. In wing patterns and shape, resembles a sympatric hairstreak Chrysoritis lycegenes (Trimen, 
1874) (possible mimicry), i.e. wings are more rounded than most Lycaena and forewing black spots are 
closer to the margin, hindwing without tails and patterned more similar to Polyommatinae than to most 
Lycaena : pale-brown with darker marginal lunules and with paler spots usually darker in the middle. 

Etymology. The name is a masculine noun in the nominative singular, formed as L[ycaen]a + [A]/f[ica] + 
on to indicate African origin of the genus and reach gender agreement with the type species name. 

Species included. The type species and Lycaena clarki Dickson, 1971, both from South Africa. 

Parent taxon. Subfamily Lycaeninae [Leach], [1815]. 
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Lycaeninae genera, subgenera and their available synonyms 


Here, we update the Appendix of Sibatani (1974) and suggest the following treatment of Lycaeninae 
grouped into nine genera. Placements of Lafron Grishin, gen. n. and Phoenicurusia Verity, 1943 are 
provisional due to the lack of both genomic data and unambiguous phenotypic evidence, and follow 
published works based on morphology and limited DNA analysis (Klots 1936; Sibatani 1974; van Dorp 
2004; de Jong and van Dorp 2006). The list is preliminary and further changes are expected in groups 
poorly covered by our genome-based phylogeny. Junior subjective synonyms are preceded by "=". 
Unavailable names are not listed. Type species are given in parenthesis with their original genus name. 

Genus Lafron Grishin, gen. n. ( Papilio orus Stoll, [1780]) 

Genus Lycaena [Fabricius], 1807 (Papilio phlaeas Linnaeus, 1760) 

Subgenus Lycaena [Fabricius], 1807 {Papilio phlaeas Linnaeus, 1760) 

Subgenus Thersamolycaena Verity, 1957 {Papilio dispar Haworth, 1802) 

Subgenus Heodes Dalman, 1816 {Papilio virgaureae Linnaeus, 1758) 

=Loweia Tutt, 1906 {Papilio dorilis Hufnagel, 1766) 

=Thersamonia Verity, 1919 {Papilio thersamon Esper, 1784) 

=Palaeochrysophanus Verity, 1943 {Papilio hippothoe Linnaeus, 1760) 

=Alciphronia Kogak, 1992 {Papilio alciphron Rottemburg, 1775) 

=Mirzakhania Kogak, 1996 {Chrysophanus kasyapa F. Moore, 1865) 

Genus Helleia Verity, 1943 ( Papilio helle Denis & Schiffermuller, 1775) 

Genus Tharsalea Scudder, 1876 ( Polyommatus arota Boisduval, 1852) 

Subgenus Epidemia Scudder, 1876 (Polyommatus epixanthe Boisduval & Le Conte, [1835]) 

=Hyllolycaena L. Miller & F. Brown, 1979 (Papilio hyllus Cramer, 1775) 
llellolycaena Ko?ak, 1983 (=Polyommatus thoe Guerin-Meneville, [1832], which is Papilio hyllus Cramer, 1775) 

Subgenus Chalceria Scudder, 1876 ( Clvysophanus rubidus Behr, 1866) 

=Gaeides Scudder, 1876 (Chrysophanus dione Scudder, 1868) 

Subgenus Tharsalea Scudder, 1876 (Polyommatus arota Boisduval, 1852) 

Subgenus Hermelycaena L. Miller & F. Brown, 1979 ( Chrysophanus hermes W. H. Edwards, 1870) 

Subgenus Hyrcanana Bethune-Baker, 1914 (Polyommatus caspius Lederer, 1870) 

~Sarlhusia Verity, 1943 (Polyommatus sarthus Staudinger, 1866) 

Subgenus Phoenicurusia Verity, 1943 (Polyommatus phoenicurus var. margelanica Staudinger, 1881) 

=Athamanthia Zhdanko, 1983; (Polyommatus athamantis Eversmann, 1854) 

Genus lophanus Draudt, 1920 (Chrysophanus (?)pyrrhias Godman & Salvin, 1887) 

Genus Heliopltonis Geyer, [1832] (=H. belenus Geyer, [1832], which is Polyommatus epicles Godart, [1824]) 
Subgenus Heliophorus Geyer, [1832] ( =H.. belenus Geyer, [1832], which is Polyommatus epicles Godart, [1824]) 

=Ilerda E. Doubleday, 1847 {Polyommatus epicles Godart, [ 1824]) 

=Kulua Zhdanko, 1995 {Polyommatus tamu Kollar, 1844) 

Subgenus Nesa Zhdanko, 1995 {Polyommatus sena Kollar, 1844) [not a homonym! Nesa Leach, 1818 is a misspelling] 
Genus Apangea Zhdanko, 1995 {Chrysophanus pang Obertluir. 1886) 

Genus Melanolycaena Sibatani, 1974 (Melanolycaena altimontana Sibatani, 1974) 

Genus Boldenaria Zhdanko, 1995 {Lycaena boldenarum White, 1862) 


Habrodais Scudder, 1876, Favonius Sibatani & Ito, 1942, Neozephyrus Sibatani & Ito, 
1942, Quercusia Verity, 1943, Chrysozephyrus Shirozu & Yamamoto, 1956, and 
Sibataniozephyrus Inomata, 1986 are junior subjective synonyms 

of Hypaurotis Scudder, 1876 

Inspecting genomic phylogenetic trees of US butterfly species (Zhang et al. 2019d), we noticed a close 
relationship between the only two New World genera from the tribe Theclini Swainson, 1830: Hypaurotis 
Scudder, 1876 (type species Thecla crysalus W. H. Edwards, 1873) and Habrodais Scudder, 1876 (type 
species Thecla grunus Boisduval, 1852). Despite dissimilar wing patterns and colors, divergence between 
the type species of these genera is indeed comparable to that of congeners (Fig. 6 top) and is even lower 
than the divergence in some compact genera, such as Chlorostrymon Clench, 1961, Ministrymon Clench, 
1961 or Electrostrymon Clench, 1961 (Fig. 6 top), and particularly in more diverse genera such as 
Strymon Hiibner, 1818, Callophiys Billberg, 1820, or Satyrium Scudder, 1876 (Zhang et al. 2019d). COI 
barcodes of Hy. crysalus and Ha. grunus are only 4.3% (28 bp) different: divergence similar to that at 
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times reported for different individuals of 
the same species (Zakharov et al. 2009; 

Kodandaramaiah et al. 2013), strongly 
suggesting that these two species are 
congeneric. The combination Hypaurotis 
grunus has been used previously in 
publications (Garth 1934), however, we 
failed to find the combination Habrodais 
ciysalus published. Also, //. crysalus has 
a slightly broader distribution and is a 
more familiar butterfly. Therefore, among 
the two names published in the same work 
(Scudder 1876), we use Hypaurotis as 
valid, and Habrodais as its junior subjective synonym, resulting in Hypaurotis grunus (Boisduval, 1852) 
and Hypaurotis poodiae J. Brown & Faulkner, 1982, revised and new combinations. 

Next, when a genomic dataset of an Old World species Favonius quercus (Linnaeus, 1758) was 
included in the Theclinae tree, it clustered closely with Hypaurotis (Fig. 6 top). The COI barcodes of H. 
ciysalus and F. quercus differ by 4.7% (31 bp). Because we lack genomic data for other genus-group 
names from the Thecla section of Eliot (Eliot 1973), we downloaded available COI barcode data from 
GenBank (Sayers et al. 2020). We find that while for most species pairs, e. g., H. ciysalus and Favonius 
orientalis (Murray, 1875), the barcodes are more similar (6%), for others, e. g., H. ciysalus and Thecla 
betulae (Linnaeus, 1758) (the type species of Thecla Fabricius, 1807) the difference is larger (8.2%). The 
COI barcode distance dendrogram computed using BioNJ (Gascuel 1997) as implemented by the 
phylogeny.fr server (Dereeper et al. 2008) reveals close clustering of species we propose to place in 
Hypaurotis (Fig. 6 bottom, red) and separation of T. betulae (green) from this cluster. In the absence of 
genomic data, this COI barcode analysis of their type species suggest that in addition to Habrodais 
Scudder, 1876 and Quercusia Verity, 1943 (type species Papilio quercus Linnaeus, 1758), the following 
four genus-group names Sibataniozephyrus Inomata, 1986 (type species Zephyrus fujisanus Matsumura, 
1910), Neozephyrus Sibatani & Ito, 1942 (type species given as Thecla taxila Bremer, 1861, which was, 
however, a misidentified Dipsas japonica Murray, 1875; according to the Art 70.3.2. of the ICZN Code 
the actual taxonomic identity of this species is chosen and japonica is fixed as type species), 
Chiysozephyrus Shirozu & Yamamoto, 1956 (type species Thecla smaragdina Bremer, 1861) and 
Favonius Sibatani & Ito, 1942 (type species Dipsas orientalis Murray, 1875) are junior subjective 
synonyms of Hypaurotis Scudder, 1876, which is a genus distinct from Thecla Fabricius, 1807. In accord 
with genetic similarities, all these species are similar in appearance (Eliot 1973). We expect that future 
studies will reveal additional synonyms and possibly a subgeneric structure of Hypaurotis. 

Finally, even from a practical standpoint of American butterfly knowledge, it seems more instructive 
to treat H. crysalus and II. grunus as congeneric emphasizing on their close kinship (despite apparent 
phenotypic dissimilarity), instead of placing them in two monotypic (or nearly monotypic) genera that 
accentuate their tenuous (but superficial) uniqueness. Finding their close relatives in the Old World places 
Hypaurotis among other Holarctic Theclinae genera such as Callophiys and Satyrium and emphasizes 
somewhat unusual but recurrent pattern revealing the connection between the Old and the New Worlds. 


.Hypaurotis crysalus|17ll7A05|USA:CO,Douglas Co.11993 

-i Habrodais grunus|PA0466|USA:CA,Plumas Co.|2017 
- Favonius (=Quercusia) quercus|20035A04|France:Var|2018 
i—— Chlorostrymon simaethis sarital5161|USA:TX,Starr Co.|2015 
“LJ-,Chlorostrymon maesites|18033Fll|USA:FL,Monroe Co,|2014 
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Chlorostrymon telea 110605 |Mexico:Tamaulipas| 1975 
Ministrymon azia|17119F07|USA:FL,Broward Co.,Coral Springs|1989 
Ministrymon janevicroyl 17119F06I Mexico :Sonora| 1995 
Ministrymon leda110610|USA:TX,Culberson Co,|1978 
Ministrymon dytie|4021|USA:TX,Cameron Co.|2015 
i— Electrostrymon angelia|15115G04|USA:FL,Miami-Dade Co.|1982 
lr 1 , Electrostrymon hugon|15115F02|Belize:Cayol2003 

Electrostrymon iovall5ll5F07IPanama.Canal Zonell979 


Hypaurotis crysalus|17117A05|USA:CO,Douglas Co. 

- Chry: 


Chrysoz.ephyrus s maraadinus [GU372536 
Habrodais grunus[PA0466JUSA:CA,Plum 
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£195528 
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*£195514 


Fig. 6. Hypaurotis , Habrodais and others. Nuclear genome tree above, 
COI barcode dendrogram below. Specimens without locality given are 
from GenBank and their accession numbers are indicated. 


Family Riodinidae Grote, 1895 

Plesioarida Trujano & Garda, 2018 is a 
junior subjective synonym of Roeberella Strand, 1932 

Plesioarida Trujano & Garcia, 2018 (type species Apodemia walkeri Godman & Salvin, 1886) was described 
as a genus (Trujano-Ortega et al. 2018) and was treated as a subgenus of Apodemia C. Felder & R. Felder, 
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1865 (type species Lemonias 
mormo C. Felder & R. Felder, 
1859) by Zhang et al. (2019e) to 
present a more internally consistent 
classification of the tribe Emesidini 
Seraphim, Freitas & Kaminski, 
2018. Continuing with the genomic 
sequencing of Riodinidae, we were 
surprised to find that a syntype of 
the type species of the genus 
Roeberella Strand, 1932, i.e., 

Lemonias calvus Staudinger, 1887 
(type locality Peru: Chanchamayo), 
was in the same clade with A. 
walkeri (type locality Mexico: 
Guerrero), rendering Plesioarida 
paraphyletic in all three trees that 
we routinely construct (Fig. 7). 
Statistical support for the 
placement of Roeberella calvus 
within the subgenus Plesioarida as 
sister to both A. walkeri and 
Apodemia hepburni (Godman & 
Salvin, 1886) (type locality 


nuclear genome 



Roeberella calvusl 18054E06|ST|M|Peru < Chanchamayo|prior to 1887|ZMHB 
Apodemia walkeri14514|M|USA:TX,Hidalao Co.|2015 
Apodemia hepburni hepburnil 17066E05[Mexico:Sonora|2003 
Apodemia palmerii arizonal 10226|USA:AZ,Pima Co.|2017 
Apodemia murphyil 17066D06|M|Mexico:Baja California Sur|2004 
Apodemia hypoglauca hypoglauca|18048Cll|M|Costa Rica:Guanacaste| 1973 
Apodemia ares j 17114H011 Ml USA: NM,Hidalgo Co.,Pine Canyon|1991 
Apodemia arnacisl 18045Fll|M|Mexico:Oaxaca| 1990 

Apodemia zela cleisj8217[F|USA:AZ,Santa Cruz Co.,Coronado N.F.|2017 


Apodemia chisosensisl 154|F|USA:TX,Brewster Co.,Big Bend National Park|2004 
Apodemia nais|PA0145|USA:C0,Gilpin Co.,Golden Gate Canyon S. P.|2016 

— Apodemia multiplaga|7635|F|USA:TX,Cameron Co.,Brownsville| 1974 
Apodemia virgulti virgultil 15105E04|NTlM|USA:CA,Los Angeles Co.|1951|CAS 
-z 02 | M [r-—--- 


Apodemia mormo mormo| 17107E02|M|USA:WA,Kittitas Co.11990 
* Apodemia mejicanus mejicanus|13383BrockC10|F|USA:AZ,Santa Cruz Co.11990 
Apodemia duryi|700|M|USA:TX,Brewster Co.,Big Bend National Park|2009 


Z chromosome 


Roeberella calvusl 18054E06|ST|M|Peru,Chanchamayo|prior to 1887|ZMHB 


Apodemia walkeri|4514lM|USA:TX.Hidalgo Co.|2015 
— : — Apodemia hepburni hepburnil 17066E05[Mexico:Sonora|2003 


j— i Apodemia palmerii arizonal 10226|USA:AZ,Pima Co.|2017 
1 — Apodemia murphyil 17066D06|M|Me.xjco:Baja California Sur|2_004 
toglai ' ' 


Apodemia hypoglauca hypoglaucal 18048Cll|M|Costa Rica:Guanacaste| 1973 




og 

Apodemia ares| 17114H01|M I USA: NM,Hidalgo Co.,Pine Cariyonl 1991 
Apodemia arnacisl 18045Fll|M|Mexico:Oaxaca| 1990 


Apodemia zela cleis|8217|F|USA:AZ,Santa Cruz Co,,Coronado N.F.|2017 


• Apodemia chisosensisl 154|F|USA:TX,Brewster Co.,Big Bend National Parl.|2004 
Apodemia nais|PA0145|USA:C0.Gilpin Co., Golden Gate Canyon S. P.|2016 

Apodemia multiplaga|7635|F|U5A:TX,Cameron Co.,Brownsville|1974 


Apodemia virgulti virgultil 15105E04|NT|M|USA:CA,Los Angeles Co.|1951|CAS 
Apodemia mormo mormo| 17107E02|M|USA:WA,Kittitas Co.11990 

Apodemia mejicanus mejicanus| 13383BrockC10|F|USA:AZ,Santa Cruz Co.11990 
Apodemia duryi|700|M|USA:TX,Brewster Co.,Big Bend National Park|2009 


mitochondrial genome 


Roeberella cajvusl 18054E06|ST|M|Peru,Chanchamayo|prior to 1887|ZMHB 



Apodemia walkeri14514|M|USA:TX,Hidalgo Co.|2015 
Apodemia hepburni hepburnil 17066E05|Mexico:Sonora|2003 

- Apodemia palmerii arizonafl0226|USA:AZ,Pima Co.|2017 

- Apodemia murphyil 17066D06|M|Mexico:Baja California Sur|2004 
- Apodemia hypoglauca hypoglaucal 18048Cll|M|Costa Rica:Guanacaste|1973 

Apodemia ares| 17114H01|M|USA:NM,Hidalgo Co..Pine Canyon|1991 
Apodemia arnacisl 18045F111 M|Mexico:Oaxaca) 1990 

— Apodemia zela cleis|8217fF|USA:AZ,Santa Cruz Co.,Coronado N.F.|2017 
Apodemia chisosensisl 154|F|USA:TX,Brewster Co., Big Bend National Park|2004 
Apodemia nais|PA0145|USA:C0,Gilpin Co., Golden Gate Canyon S. P.|2016 

- Apodemia multiplaga|7635|F|USA:TX,Cameron Co.,Brownsville| 1974 

3 Apodemia virgulti virgultil 15105E04|NT|M|USA:CA,Los Angeles Co.|1951|CAS 
Apodemia mormo mormoll7107E02|M|USA:WA,Kittitas Co.1 1990 
Apodemia mejicanus mejicanusl 13383BrockC10|F|USA:AZ,Santa Cruz Co.1 1990 


Apodemia duryi|700|M|USA:TX,Brewster Co.,Big Bend National Park|2009 

Fig. 7. Subgenus Roeberella (red), and other Apodemia (blue). 




A. h. hypoglauca 


A. (Roeberella, =Plesioarida ) calvus 


syntype; 
type species 


Fig. 8. Previously unnoticed but apparent phenotypic similarity between Roeberella calvus (c, syntype and the type species 
of Roeberella, Peru: Chanchamayo, NVG-18054E06 in ZMHB, here placed in Apodemia), and representative species of the 
subgenus Plesioarida (here synonymized with Roeberella)'. Apodemia walkeri (d, the type species of Plesioarida, Mexico: 
Morelos, NVG-18048D04 in USNM) and Apodemia hypoglauca (a, Costa Rica, NVG-18048C11; b, Mexico: Morelos, 
NVG-18048C10 both in USNM). Dorsal above, ventral below. This similarity was uncovered by genomic analysis. 


NVG-18048C11 


NVG-18048C10 


NVG-18048D04 
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Mexico: Chihuahua) is strong in the protein coding regions of the entire nuclear genome, Z-chromosome 
and mitochondrial genome: 100% of all segment trees contained the clade of these three species. Such a 
result was unexpected, because the type species of Roeberella (South American) was not previously 
compared with Apodemia (North American), and we suspected a possibility of error or contamination. 
However, COI barcodes of R. calvus and A. walkeri differ by only 5.2% (34 base pairs), a difference 
smaller than that between A walkeri and A. mormo (8.5%, 56 bp). Moreover, R. calvus is phenotypically 
similar to Plesioarida species (Fig. 8) and shares a falcate forewing with Apodemia hypoglauca (Godman & 
Salvin, 1878) (type locality "Mexico"). The similarities are most prominent in the ventral wing pattern (Fig. 
8bc bottom). Interestingly, the wing shape of R. calvus appears more like that of female Plesioarida (Fig. 
8b). Dorsally, R. calvus reminds us of an aberrant Plesioarida. Therefore, due to genetic and phenotypic 
similarities, we conclude that Plesioarida syn. n. is a junior subjective synonym of Roeberella, and 
consequently we place Roeberella as a subgenus (new status) of Apodemia. 

We think this find is particularly interesting for several reasons. First, it extends the range of 
Apodemia, previously not recorded from South America, to Peru, with all evolutionary and 
biogeographical implications of this fact. Second, it underscores the importance of a more comprehensive 
phylogenetic analysis before proposing new genus-group names to avoid creation of unnecessary 
synonyms (Trujano-Ortega et al. 2018; Trujano-Ortega et al. 2020). Third, it reiterates the power of 
genomic approaches and the value of the type concept, both in species and genus-group names. 
Sequencing of the syntype of the type species of the genus-group name Roeberella solidifies our 
conclusions, eliminating a possibility of misidentification or incorrect inference from a non-type species. 


Lemonias lencates Hewitson, 1875 currently placed in Roeberella Strand, 1932 

belongs to Pachythone H. Bates, 1868 

Genomic sequencing of a specimen of Roeberella lencates (Hewitson, 1875) (type locality not given, 
species recorded from Costa Rica to Brazil) and comparisons among the type species of available genus- 
group names reveal that R. lencates (Fig. 9, red) is not monophyletic with the type species of Roeberella 
Strand, 1932 (Lemonias calvus Staudinger, 1887, type locality Peru: Chanchamayo, a syntype sequenced, 
Fig. 9, magenta), but instead belongs to 
Pachythone H. Bates, 1868 (type species 
Pachythone erebia Bates, 1868) (Fig. 9, 
blue). The tree includes the type species of 
all available genus-group names in 
Emesidini placed among other Riodinidae 
(Fig. 9, taxa currently placed in Apodemia 
are in green font). Therefore, we transfer R. 
lencates and a number of its close allies as 
suggested by their phenotypic similarities 
from Roeberella to Pachythone, forming the 
following new combinations: Pachythone 
lencates (Hewitson, 1875), Pachythone flocculus (Brevignon & Gallard, 1993), Pachythone floccus 
(Brevignon, 2013), Pachythone heberti (P. Jauffret & J. Jauffret, 2007) and Pachythone marajoara (P. 
Jauffret & J. Jauffret, 2007). Furthermore, due to close clustering in the genomic tree (Fig. 9) comparable 
to that within other genera (e.g., Periplacis Geyer, 1837 in addition to Emesis Fabricius, 1807 and 
Apodemia) and phenotypic similarities in wing patterns (e.g., compare with Pachythone strati (Kaye, 
1925) and Pachythone rubigo (H. Bates, 1868)), we suggest that Pseudonymphidia Callaghan, 1985 (type 
species Emesis clearista Butler, 1871. Fig. 9, cyan) is a subgenus (new status) of Pachythone (Fig. 9, the 
clade labeled with the name in red font highlighted in yellow). 


Pachythone 


Pachythone erebia {Pachythone}|19129C01|Ecuador|2015 
Pachythone strati (=corculum) {Pixus}| 19031A09|Panama| 1981 
- Pachythone velazquezi {Lamphiotes}! 19031A10|Mexico:Gro| 1977 

Roeberella lencatesl 19029H05|Panama| 1970 
Pseudonymphidia clearista {Pseudonympnidia}| 19031A07|Panama| 1971 

- Stalachtis phaedusa {Stalachtis}|19034Gll|Guyana|2001 

- Stalachtis phlegia {Nerias}|19034G07|Brazil:RJ|1994 

Periplacis glaucoma {Peripfacis}|19033G02|Brazil:SC|1990 
Periplacis menander {Menander}|18072E04|Costa Rica|2013 
Emesis cereus {Emesis}! 18045C06|Brazil:RO| 1993 
Emesis fatimella {Nelone}ll8044G08lVenezuelal 1984 
Emesis tenedia {Tenedia}|18044H06|Costa Rica|2010 
Emesis brimo brimo {Brimia}| 18045E01|Costa Rica|2008 
Emesis mandana {Mandania}|18044D09|Guyana|1999 
Emesis lucinda {Aphacitis}| 18044B02|Guyana|2000 
Emesis poeas {Poeasia}|18044G02|Mexico:Oax|1980 
— Curvie emesia {Curvie}|5245|USA:TX|2015 
Roeberella calvus {Roeberella}| 18054E06|ST|Peru|before 1887 
Apodemia walkeri {Plesioarida>|4514|USA:TX|2015 
Apodemia mormo mormo {Apodemia}| 17107E02|USA:WA| 1990 
- Apodemia nais {Neoapodemia}| 17066F02|USA:AZ| 1993 
Riodina lysippus {Riodina}|19026H04|Ecuador|1988 

Fig. 9. Emesidini and Pachythone among others. Genus-group names 
each species is the type of are given in brackets. 
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Family Nymphalidae Rafinesque, 1815 

Speyeria Scudder, 1872 is a subgenus of Argynnis Fabricius, 1807 

A close relationship between the New World genus Speyeria Scudder, 1872 (type species Papilio idalia 
Drury, 1773) and the Old World genus Argynnis Fabricius, 1807 (type species Papilio paphia Linnaeus 
1758) has been suggested (Simonsen 2006; Simonsen et al. 2006). In these studies, Simonsen proposed to 
treat Speyeria as a subgenus of 
Argynnis. Although this suggestion 
has been followed in a number of 
works (Wells et al. 2011; Scott and 
Fisher 2014), it has not been 
universally accepted (Pelham 2008; 

De Moya et al. 2017; Pelham 
2020), likely due not to scientific 
but historical reasons. Several 
generations of American naturalists 
were raised being accustomed to 
the name Speyeria and are less 
familiar with the name Argynnis, 
thus being resistant to abandoning 
Speyeria as a genus name. Our ge¬ 
nomic studies also support the view 
that Speyeria should be considered a subgenus within Argynnis. The two groups are very close to each 
other genetically (Fig. 10). The COI barcode difference between the type species of Speyeria and 
Argynnis is 8.2%, and the estimated time of divergence is 9.1 Mya according to De Moya et al. (2017) 
(Brenthis Hubner, [1819], which is sister to the clade that includes Argynnis and Speyeria diverged from 
them 11.1 Mya), but likely about 7.5 Mya according to Chazot et al. (2019), who show that Brenthis 
diverged from Argynnis at 9.2 Mya (9.1 *9.2/11.1—7.5). This divergence is nearly the same as (or less 
than) between the two species of Euptoieta E. Doubleday, 1848 (Fig. 10): COI barcode difference 8.8%, 
estimated time of divergence about 8 Mya according to Zhang et al. (2019d). Moreover, Argynnis is also 
quite close to Boloria Moore, 1900 (type species Papilio pales [Denis & Schiffermuller], 1775) 
(estimated divergence 14.6 Mya (Chazot et al. 2019)), and they together form a more prominent group in 
the phylogenetic tree than either of them does separately (Fig. 10). Thus, it is even conceivable to take the 
next step and treat Boloria sensu lato as a subgenus of Argynnis. Currently we refrain from their 
unification, because the genetic distance between Argynnis and Boloria is still within the limits possible 
for distinct genera, and the pronounced phenotypic distinction between these two genera exists making 
their visual recognition straightforward. However, butterfly classification would be more inconsistent if 
Speyeria stays a genus distinct from Argynnis. Therefore, we agree with Simonsen et al. (2006) and place 
Speyeria as a subgenus within Argynnis. 



11Germany:Bavaria,near Graswang|2017 
?ll9077C06|Japan,Kyoto,Mount Hiei| 1947 
/C04|.Korea,SpbaeLsan| 1985 
.n]Switzerland,Furl a Pass]2017 


y^Mi’gynnis paphia|PAOE2 
' Argynnis sag[ana lianel- 
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Fig. 10. Argynnis (blue) and its subgenus Speyeria (red). 


Argynnis irene Boisduval, 1869 and Argynnis nausicaa W. H. Edwards, 1874 are 
species distinct from Argynnis hesperis W. H. Edwards, 1864 

Argynnis atlantis W. H. Edwards, 1862 (type locality USA: New York, Green Co., mostly eastern in 
distribution) and Argynnis hesperis W. H. Edwards, 1864 (type locality USA: Colorado, Jefferson Co., 
mostly western in distribution) fonn a species complex that requires further investigation (Dunford 2009). 
We obtained whole genome shotgun sequences for nearly all its taxa considered valid by Pelham (2020). 
A number of these taxa were represented by their primary type specimens to ensure correct application of 
their names (Fig. 11, indicated as HT for holotype and LT for lectotype). Because protein-coding regions 
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of species in the subgenus Speyeria are quite conserved in their sequences, in order to provide better 
discrimination between taxa, we used all genomic sequences mapped to Argynnis (Speyeria) diana 
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Fig. 11 .Argynnis irene (red), hesperis (magenta), atlantis (blue) and nausicaa (green). 
Trees constructed from nuclear genomic regions from (a) autosomes and (b) Z chromosome. 


reference genome (Zhang et al. 2019d). Furthermore, we considered the autosomes and Z chromosome 
separately due to their distinct roles in evolution and differences in resistance to introgression (Cong et al. 
2019b). The trees constructed from concatenated coding and non-coding regions of autosomes (Fig. 11a) 
and Z chromosome (Fig. lib) revealed evolutionary complexities of the taxa included. First, as expected, 
the entire group consisting of A. hesperis (Fig. 11 green, red and magenta) and A. atlantis (Fig. 11 blue) is 
monophyletic in both trees. Second, in both trees, A. hesperis is paraphyletic with respect to A. atlantis, 
and a clade composed of south-central subspecies of A. hesperis (Fig. 11, green), including A. hesperis 
nausicaa W. H. Edwards, 1874 (type locality USA: Arizona, Graham Co.) is sister to all other taxa 
combined. Third, the position of the clade consisting of northern and eastern subspecies of A. hesperis 
that include the nominotypical subspecies (Fig. 11 magenta) is different in autosome and Z chromosome 
trees. In the autosome tree (Fig. 1 la), this clade is sister to A. atlantis, and in the Z chromosome tree (Fig. 
lib), it is sister to the clade of western subspecies of A. hesperis that include A. hesperis irene 
(Boisduval, 1869) (type locality USA: California, Sierra Co.). 

Thus, the trees reveal four groups of taxa in this complex, and applying the oldest name in each 
group, we call them atlantis, hesperis (sensu stricto), irene and nausicaa. To probe whether these groups 
are species, we used the Fst/Gmin Z chromosome tests (Cong et al. 2019b) obtaining the following 
statistics. First, traditionally treated as distinct species, atlantis and hesperis groups show differences 
consistent with their species-level distinction: 0.23/0.07, albeit marginally (Fst for distinct species is 
typically above 0.2, with 0.5 and above indicating strong differentiation, and Gmin is less than 0.1, with 
0.02 and below indicating strong isolation). Second, the differences between the atlantis and nausicaa 
groups are even more pronounced than those for atlantis and hesperis : 0.36/0.026, which is in agreement 
with the nausicaa group being sister to the clade consisting of the three other groups in both autosome 
and Z chromosome trees (Fig. 11). Third, the differences between hesperis and nausicaa are of about the 
same magnitude as for others: 0.26/0.05, indicating that A. nausicaa is a distinct species rather than a 
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group of subspecies within paraphyletic A. hesperis. Fourth, the hesperis and Irene groups did not reveal 
species-level differences in the Z chromosome: 0.13/0.13, suggesting that they may be conspecific, which 
is in agreement with their close clustering together in the Z chromosome tree (Fig. 1 lb). 

Although the irene group is not strongly different from the hesperis group in Z chromosome, it is 
placed differently in the autosome tree: as sister to both hesperis and atlantis groups, rendering the 
species A. hesperis that includes irene paraphyletic. While a species paraphyletic in a tree built from 
concatenated genomic alignments is not inconceivable due to the possibility of extensive introgression 
from some other species in a part of the species range, such a situation calls for further investigation. 
Analysis of the trees built from various segments of the nuclear genome revealed that some segments in 
the hesperis group are similar to the atlantis group, while other segments are similar to the irene group. 
Hence, we hypothesize that A. hesperis is a hybrid species of A. irene and A. atlantis, because it shares 
20% and 67% of its autosome-linked genome, and 71% and 20% of its Z-linked genome with the latter 
two species respectively, while it possesses only 0.17% of unique polymorphisms, compared to 0.5% and 
0.28% unique polymorphisms in A. irene and A. atlantis, respectively. We see that a significant fraction 
of the A. hesperis genome is shared with either A irene or A. atlantis, and the number of unique mutations 
in the A. hesperis lineage is smaller than that compared to either of its putative parental species, 
suggesting a hybrid origin of A. hesperis. Consequently, we consider A. irene (consisting of four 
westernmost subspecies presently associated with A. hesperis. Fig. 11) to be a species-level taxon, and the 
Z chromosome similarity with A. hesperis is therefore explained by the hybrid origin of A. hesperis, 
which inherited larger segments of this chromosome from A. irene. This scenario of species originating by 
hybridization is not covered by the Fst/Gmin Z chromosome test for species distinction (Cong et al. 
2019b). The COI barcodes of A. irene are closer to A. atlantis (2.5%, 17 bp difference) than to A. hesperis 
(5%, 33 bp difference), probably because A. hesperis possesses mitogenomes introgressed from A. 
nausicaa and does not reveal differences in the barcodes with the latter species. As a side note, the earlier- 
named species A. hesperis is likely to a be a hybrid species with one of the parental species being a named 
later (A. irene), illustrating that biological reality has little to do with the order species were named in. 

Finally, we find (Fig. 11) that the holotype of Speyeria hydaspe conquista dos Passos & Grey, 1945 
(type locality USA: New Mexico, Santa Fe Co., presumed to be in error), presently placed as a synonym 
of Argynnis hydaspe rhodope W. H. Edwards, 1874, clusters closely with the holotype of Argynnis 
hesperis tetonia (dos Passos & Grey, 1945) (type locality USA: Wyoming, Teton Co.) and is therefore 
placed as a synonym of tetonia, new placement. 

In summary, genomic data suggest that the atlantis-hesperis complex consists of four species: A. 
atlantis, A. hesperis, and two others with reinstated status: A. irene and A. nausicaa. The following 
subspecies are assigned to A. irene to form new combinations: Argynnis irene dodgei Gunder, 1931, 
Argynnis irene cottlei J. A. Comstock, 1925, and Argynnis irene hanseni (J. Eminel, T. Emmel & 
Mattoon, 1998). The following subspecies are assigned to A. nausicaa to form new combinations: 
Argynnis nausicaa elko (Austin, 1984), Argynnis nausicaa greyi (Moeck, 1950), Argynnis nausicaa viola 
(dos Passos & Grey, 1945), Argynnis nausicaa tetonia (dos Passos & Grey, 1945), Argynnis nausicaa 
chitone W. H. Edwards, 1879, Argynnis nausicaa schellbachi (Garth, 1949), Argynnis nausicaa electa W. 
H. Edwards, 1878, Argynnis nausicaa dorothea (Moeck, 1947), and Argynnis nausicaa capitanensis (R. 
Holland, 1988). The names for other taxa in this complex remain unchanged. 


Argynnis zerene atossa W. H. Edwards, 1890, new combination, 
is not a subspecies of Argynnis adiaste W. H. Edwards, 1864 

Described as a species Argynnis atossa W. H. Edwards, 1890 (type locality USA: California, Kern Co.), 
this likely extinct butterfly was placed as a subspecies of Argynnis adiaste W. H. Edwards, 1864 (type 
locality USA: California, Santa Cruz Co.) due to its prominent wing pattern similarities, in particular, the 
washed out ventral hindwing devoid of silvery spots and reduced black markings above, especially on its 
hindwing. To our surprise, none of the trees (concatenated protein-coding regions of the entire nuclear 
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genome, of Z chromosome and mitochondrial genome) placed atossa with adiaste (Fig. 12). Instead, the 
atossa clade originated within Argynnis zerene Boisduval, 1852 (type locality USA: California, possibly 
Plumas Co.), which is phenotypically strongly spotted, both in back above and silvery below. Argynnis 
adiaste was sister to Argynnis hvdaspe Boisduval, 1869 (type locality USA: California, possibly Sierra 
Co.) in a clade remote from A. zerene. Using Fst and Gmin statistics, atossa and A. adiaste adiaste are 
characterized by 0.48 and 0.012, respectively, indicating strong genetic isolation. However, atossa and A. 
zerene show the values Fst=0.13 and Gmin=0.15, which are within the range characteristic of 
conspecificity. For comparison, A. adiaste adiaste and A. adiaste clemencei J. A. Comstock, 1925 (type 
locality USA: California, San Luis Obispo Co.) exhibit Fst/Gmin of 0.06/0.21, showing less genetic 
differentiation from each 
other than atossa from A. 
zerene. COI barcodes of 
atossa and adiaste adiaste 
lectotypes differ by 2.7% 

(18 bp), but the A. zerene 
myrtleae (dos Passos & 

Grey, 1945) holotype and 
atossa lectotype differ by 
1.2% (8 base pairs). We 
think that the evidence 
presented here to support 
that atossa is not A. adiaste 
is strong. However, while 
currently we do not have 
data to justify it, a 
possibility that atossa is a 
species distinct from A. 
zerene exists, if A. zerene is 
later found to be a complex 
of several distinct species. 

For the lack of a better 
option, we place atossa as a 
subspecies of A. zerene, 
where it fits genetically not 
worse than a number of 
other A. zerene subspecies, 
although it differs strongly from all other A zerene subspecies in wing patterns. 

The association of atossa with zerene came as a surprise due to their phenotypic dissimilarity. The 
reasons behind the wing pattern similarities between atossa and adiaste, and the lack of such similarity 
between atossa and zerene remain unclear. We sequenced the primary types of both atossa (in the 
Carnegie Museum of Natural History, collected prior to its description in 1890) and tejonica J. A. 
Comstock, 1925 (in the Los Angeles County Museum, considered a subjective junior synonym of atossa, 
has spots weakly silvered, collected in 1923 in Los Angeles County, California) in addition to a specimen 
from Colorado State University collection (collected in 1922 in Sierra Madre Mountains, California). 
These specimens were collected in different localities and different years, and handled differently 
throughout the years, but they cluster together in all trees, forming a distinct clade within zerene. It does 
not seem likely that legs of yet unsequenced population of zerene were glued to all these 3 specimens, or 
similar contaminant affected these samples processed in our lab on different days (and for some of them, 
years) and sequenced in different lanes and batches. Therefore, we conclude that Argynnis zerene atossa 
W. H. Edwards, 1890, new placement, represents an unusual example of phenotypic convergence that 
hindered its taxonomic placement, now revealed through genomic sequencing. 
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- Speyeria idalia occidental^ 1/115A04|M|.USA:CO,Kit Carson Co. A miles east of Flagler|1987 

Speyeria diana|9190 MlUSAiAR.Polk Co.,Ouachita National Forest]2017 
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Speyeria cybelecybele^lPlLFjUSAiAR.Scott Co..Ouachita National Forest|2017 . 
Speyeria aphrodite etnne1PAO220 IUSA: CO,Larimer Co.,2 miles NNE Virginia Dale|2016 
Speyeria nesperis nausiqaaj 11628LM|,USA:AZ.Greenlee Co.12018 
peyeria atlantis atlantisri9065C03 Canaaa:Ontario|2013 

• - 1 - ---ark!2017 


Speyeria.mormoma e.urynomelOSSOjUSA:WY.Park Co.,Yellowstone National Parkl2( 

- Speyeria edwardsii|PA0218|USA:C0,Larimer Co.,2 miles NNE Virginia Dale[2016 

- Speyeria hydaspe minpr[01506CQ5lHT[M|Canada:Britisn Columbia,Lillooet, Mount McLeal 1926 

Speyeria hydaspe tenama|l9061A07lHT MIUSAiCA,T ehama Co.,Covelo Roach at Government |1965 
Speyeria hydaspe rhodope|9431|USA:WY,Park Co..Yellowstone National Park|2017 

-u. j -l SA: wY y P ar k Cp..Yellowstone 1 ■ — 

tilla. Co.,Umatilla 


Speyeria hydaspe rl 


'6301'fi [USA : 6 k ; r Um atill: 


National Park|2,017 
a National Forest|2017 
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— Speyeria hydaspe tehama]PA0418lySA:CA,Mendocino Co..Mendocino Pass|2i 

Speyeria hydaspe nyaaspell8669A12[LT|M|USA:CA,Sierra Co.,feold Lake|priorto :_ 

■ Speyeria hydaspe.hy.daspe (=purpurascens)J19044B01|LT[M|USA:CA,Siskiyou CoJprior to 1 
Speyeria hydaspe viridicormsl 19061A061HTJM USA:CA,Kern Co.,Greenhorn Mts.|1923 
*ria hydaspe davisillOOeiAOSLHTlMJklSAiC^Humboldt Co.,McCleJIan Mountain^* 


Speyeria'hydaspe davisi|1906iACL „ 

• peyeria adiaste adiaste 119051 


_,_..Janta Cruz Co.,Saratoga Gapll_ _ 

Speyeria adiaste adiaste[19051Hl?|F|U5A:CA,Santa Cruz Co..Saratoga Gap|1959 , 
Speyeria adiaste aaiasteJ19062F03TM| USA: CA.Santa Cruz Co.,NlW of Boulder CreeL|1960 
Speyeria adiaste adiastefl9051Hll|F|USA:CA,Santa Cruz Co.,Saratoga Gap|1959 
Speyeria ^diaste ^diast^ellOOeSBl^USAiCA^Santa Cruz Cp.1198'' 


j>eatb|1987 


Speyeria adiaste adiaste] 19065C0l[USA: CA.Santa Cruz Co!| 1982 
Speyeria adiaste adiaste119051H10|,M|USA:CA,Santa Cruz Co..Saratoga Gap|1959 , 

Speyeria adiaste adiaste] 19062F02]M|USA:CA.Santa Cruz Co..Santa Cruz Mountains 1990 
peyeria adiaste adiaste119D55C1 ljLTJMJUSA:CA.Santa Cruz Co..Santa Cruz|priorto 1864 
Speyeria adiaste demenceiri9061Fi0[MjUSA:CA.Monterey Co..Chews Ridge|1990 

Speyeria adiaste.clemenceil 190.6 lC05]HT| Ml USA,: CA.San Luis Obispo. Co. .Atascadero 11922 
— Speyeria adiaste clemenceiTl906iC06^TfFjUSA:CA,San Luis Obispo Co//Atascadero fi.922 

*r~""- Speyeria .adiaste clemenceil 19065G06jFlUSA:CA,Kern/Los Angeles Cos.,Tejpn Rangel 1930 

1 - Speyeria adiaste clemenceiri9065G05|MJUSA:CA,Kern/Los Angeles Cos.,Tejon Range] 1930 

.j—S3- Speyeria adiaste clemenceij 1906lFll|M|USA:CA,Monterey Co^Chews Ridge|1984 

-P—- S P e v eri a a d i a ste cl e m e n ce 1 11 9 9 6 5 C 0 2] U s A: C A, Al a in e d a C o . 119 73 

l_i— ss- Speyeria adiaste clemenceil 19062F05[M USA:CA.Monterey Co,.Chews Ridge[1977 

- Speyeria adiaste clemenceij 19062Fl?lM]USA:CA,Monterey Co...Chews Ridge 19.84 

Speyeria adiaste clemenceil 19062F09jMrusA:CA, Monterey Co.,2 mi N Chews Ridge|1994 
Speyeria adiaste clemenceiI.PAO30JUSA:CA.Monterey Co..Chews Ridge|2.0,16 , 

Speyeria adiaste cTemencei 19062F10 F]UbA:CA,Monterey Co..Chews Ridge|1990 


:lemencei]i9062Fll| Ml lisa: CA.Monterey Co..Chews Ridge 1990 

■‘ :UT,Davis Co,.Wasatch National Forest 2017 

2|FlUSA:UT,Davis Co..Wasatch National Forestl2017 

__1115A061MUJSA;NV.Clarl: Co..Spring Mountains|l965 

Speyeria carolae]l9043H01|HT|M|USA:NV.Clark Co..Charleston Park|priorto 1942 
Speyeria.carolael 18028A02jM|ySA:NV / Clark Co.,Spring Mountains! 1967 
- Speyeria carofaeri8028A03|FlUSA:NV,Clark Co..Spring Mountains] 1967 
Speyeria coronis hennei11,9044A03|HT|M[USA:CA.Ventura Co.,Mt. Pinos|prior to 193.4 , 

Speyeria coronis simaethal 19044B05]HT|M|USA:WA,Okanogan Co..Cascade Mountains|before 1945 



enemalcolmi] 19061C 121 HT|MfuSA:CA, Mono Co.,Mammoth| 1918 , 

Speyeria zerene malcolmi|PA0467|USA:CA,Mono Co.,microwave above Conway Summ|2017 
• Speyeria zerene zereneIPAO460JUSA:CA,Plumas Co..Mohawk-Chapman Road|2017 
Speyeria zerene z.ereneJPAO90|U i SA:CA / Sierra Co.,1 mi E of Calpine|2016 i 
Speyeria zerene behrensnlPA0419|USA:CA.Mendocino Co..Mendocino Pass|2017 

— Speyeria adiaste atossaLl8037H09jMjUSA:CA,Santa Barbara or.Ventura Cos. .Sierra M|1922 
Speyeria adiaste atossa]l9055D01[LTjFluSA:CA,Kern Co.,Tenacnapi| 1889 
peyeria adiaste atossa (=tejonica)Ll9061A02[HT|FlUSA:CA,Los Angeles Co.,Tejon region|1923 
Spevena zerene glonosal 19043 Hi2|HTTM|USa:OR,J osephine Co.,Illinois River vallevlpnor to 1957 
— Spey.eria zerene zerene (=conchyliatus)|19061C07|HT|M|USA:CA,Siskiyou Co,,Mt. Shasta|1921 
- t, --e|l9044A10jHTrM|USA:CA,San Mateo Co.,^an Mateolp---"* 


_ _zen 

Speyeria zerene myrtleae 
Apeyeria zerene puntareyes 
ipeyeria zerene myrtleael" 1 ' 




«,San Mateo Co.J>an Mateo]prior to 1945 
SA:CA,Marin Co.Point Reves Peninsula|1964 
■iiwjM.^M.Marin Co..Point Keyes 1972 
HTjM|USA:CA,Sonoma Co.11991 


Speyeria zerene sonomensis|19061ClOjm «-u.| 

Speyeria zerene gunderij 19061C08lHTlML.USA:CA,Mqaoc Co.j 1924 
Speyeria zerene gunderi (=cynna)|19o43H07 HUM]USA:NV,Elko C 
Speyeria zerene picta (=garretti)|19043Hll[HT M[USA:Canada,r 

Speyeria zerene pIatinaT=pfbutsi)J19044All|HTjM|USA:UT,Utal. l^ici iuc.u.c x: 

Speyeria zerene sinope [19044BQ6LHT|M|USA;CO,Larimer Co.,Rocky Mountain NPIbefore 
Speyeria zerene picta|01506B10]HT|M|Canada:British Columbia,Aspen Grove]1922 
Speye.ria zerene sinope 164671 Ml USA: CO.Gr.and Co.,4,3 ajr mi S.E Hot Sulphur SpJ2(Jl6 
Speyeria zerene picta|9626|M|USA:OR,UmatiTla Co..Umatilla National Forest|201/ 


...i Co..Humboldt NFlbefore 1945 
ida,British Columbia|before 1924 
Utah Co..Mt. Loafer before 1933 
.Rocky Mou ' ‘ 

Dia^\spen ( 


Fig. 12. Argynnis zerene atossa (red) is not A. adiaste (cyan), but A. zerene (green). 

































































































Argynnis coronis carolae dos Passos & Grey, 1942 is a subspecies-level taxon 


Described as a subspecies of Argynnis (Speyeria) coronis Behr, 1864 (type locality USA: California, 
possibly Santa Clara Co.), Argynnis coronis carolae dos Passos & Grey, 1942 (type locality USA: 
Nevada, Clark Co.) was elevated to species by Emmel and Austin (Emmel and Austin 1998) due to a 
number of its unique features, including the isolated locality. Our genomic trees revealed that carolae 
forms a compact clade consistent with its isolation, however, not separately, but deep within A. coronis, 
sister to Argynnis coronis hennei Gunder, 1934 (type locality USA: California, Ventura Co.) and near the 
nominal A. coronis (Fig. 12 magenta inside blue clade). This placement in the tree is geographically 
meaningful, but makes it difficult to accept the species status of carolae because it renders coronis non- 
monophyletic. Furthermore, Fst/Gmin for carolae vs. coronis are 0.23/0.11. Fst is somewhat elevated due 
to genetic closeness of individuals within apparently strongly inbred carolae population, but the gene 
flow between carolae and coronis is more than two times higher than the 0.05 threshold characteristic of 
different species. The COI barcodes of the A. carolae holotype and A. coronis differ by -0.5% (3 bp). 
Therefore, we reinstate this taxon as a subspecies: Argynnis coronis carolae dos Passos & Grey, 1942. 


Aglais Dalman, 1816 and Polygonia Hiibner, [1819] 
are subgenera of Nymphalis Kluk, 1780 

In agreement with others (Opler and Malikul 1992; Layberry et al. 1998; Savela 2020), we propose that 
Aglais Dalman, 1816 (type species Papilio urticae Linnaeus, 1758) and Polygonia Hiibner, [1819] (type 
species Papilio c-aureum Linnaeus, 1758) are better treated as subgenera of Nymphalis Kluk, 1780 (type 
species Papilio polychloros Linnaeus, 1758) rather than as distinct genera. These three distinct 
phylogenetic groups are close to each other genetically with genetic distances between them of the same 
magnitude as those for taxa considered congeneric in closely related lineages, such as Vanessa 
[Fabricius], 1807 (type species Papilio atalanta Linnaeus, 1758) (Fig. 13). The times of divergence 
between Nymphalis and Polygonia 
and between Nymphalis and Aglais 
have been estimated as -7 and -11 
Mya, respectively (Chazot et al. 

2019). COI barcode difference 
between N. polychloros and P. c- 
aureum is 6.7% (44 bp), and 
between N polychloros and A. 
urticae is 7.9% (52 bp). This 
divergence is comparable to that 
between Vanessa annabella (W. D. 

Field, 1971) and Vanessa atalanta (Linnaeus, 1758) at 6.2% (41 bp), but smaller than the divergence 
between V. atalanta and N polychloros of 9.7% (64 bp). While it is not possible establish a meaningful 
COI cutoff for the genus-level divergence, these numbers comparatively indicate genetic similarity of 
these butterflies, and they also loosely correlate with their divergence times. Inspection of the genomic 
tree (Fig. 13) reveals that the most prominent internal branches (the longest) are indeed those that support 
Vanessa and its sister clade consisting of Nymphalis, Polygonia and Aglais. The clade of the latter three 
taxa is compact (Fig. 13), prominent, and genetic divergence within it agrees with the expected 
divergence within a genus. Therefore, unification of the three genera under a single genus {Nymphalis) 
would be more consistent with how other genera (e.g. Vanessa) are classified. These showy butterflies are 
quite diverse in their wing patterns and attracted significant attention, which is likely responsible for their 
oversplit classification. Also, we confirm the expected sister relationship between Nymphalis antiopa 


J ais iolPAOE20]Germany:Bavaria,near Grasswang|2017 
ais uracaeTPAOE19 Germany:Bavaria,near Grasswang 2017 
Aglais milberti subpallida|9460[M|USA:WY,Park Co.,Yellowstone National Park|2017 
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■ Polvaonia interrogationisl6030IF]USA:TX.Dallas, Co,Dallas, Norbuck Pari |2016 
^ Polygonia commari71 17D09IH USA^ 


. ___iA:tX,Polk Co.119? 1 

Polygonia satyrus satyrus I PA042|USA:CA,Monterey Conoid Coast Road off SH1|2016 
Polygonia prognel 17115B04|M[USA:SD.Roberts Co.11981 
_ • Polygonia ore.as silenus1171 lfeBOSIMjUSAsORfYamhiJ] Co.11983. . . 

|— If .Polygonia gracilis zephyrus|PA0246|usA:C0,Larimer Co.,2 mi N Virginia Dale|2016 
Polygonia naroldiil 17115B06JMIMexico:Hidalgo,18 miles northeast of Zimapanl 1981 
I _i-. Polygonia faunus hvlaslPAQ2/0 USA:CO,Larimer Co,..Sky Ranch Road off Pingree|2016 
| n - Polygonia c-album|PAOE08JSwit:zerland,Furka Pass|2017 
I— Nymphalis l-album watsonili7115B03|M|uSA:MT,Lake Co..Flathead Lai.e| 1976 
N p- Nymphalis i .californica|9635|USAiUT,Davis Co.,Wasatch National Forest|2017 
I—krr^-Nymphalis antiopa (=Tintnern)| 18&41B 121 HtlM|USA:NY,Schoharie Co.,Scnoharie|prior to 1857 
If-Wmphalis antiopa 189801FJ USA: NM,Otero Co ..Lincoln National Forest|2017 
Lr -Nymphalis cyanomelasllbl01E02fMjMexico:Veracruz,Las Vigas 
■- Nymphalis cyanomelas[15101E03 FlMexico:Veracruz,xalapa 

Vanessa annabellaj9bl4[MjUSA:UT,Davis Co..Wasatch National Forest|2017 
Vanessa cardui|8b41JF|USA:TX,Dallas Co,Dallas|2017 


Vanessa virginiensisi3866| 
' -*■-1117" 


---,-,__6|FjUSA:TX,D,_ _ _ 

Vanessa tameamea|1^115B02jM|USA:HI,Hawaiiri974 

Vanessa atalanta rubria|9826|F]USA:TX,Dallas Co.Dallas, White Rock Lake Par|2017 

- Hypanartia |ethe|,1.7118C08JMexico.:Oaxacall987 

- Symbrenthia brabirari9122 , F10)Taiwanri973 

Fig. 13. Aglais (red) and Polygonia (magenta) are Nymphalis (blue). 


Dallas Co.,Dallas, White Rock Creek|2015 

.i it - ' 1 
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(Linnaeus, 1758) and N. cyanomelas (E. Doubleday, [1848]) (Fig. 13). In summary, in a move towards a 
more consistent classification, we suggest treating Aglais and Polygonia as subgenera of Nymphalis. 


Coenonympha California Westwood, [1851] is species 
distinct from Coenonympha tullia (Miiller, 1764) 

In agreement with Kodandaramaiah and Wahlberg (2009), we find substantial genetic differentiation 
between European Coenonympha tullia (Muller, 1764) (type locality Denmark: Zealand Island) and North 
American Coenonympha California Westwood, [1851] (type locality USA: California, near San 
Francisco) (Warren et al. 2016). The latter has frequently been treated as a subspecies of the former 
(Pelham 2008; Pelham 2020). 

Fst/Gmin statistics for C. tullia vs. 

California are 0.22/0.02, indicating 
very limited gene exchange between 
these taxa, typical of closely related 
species rather than subspecies. This 
differentiation strongly suggests that 
C. California is indeed a species 
distinct from C. tullia. Moreover, we 
find that the northernmost American 
taxa from this complex form a clade 
that is sister to C. tullia (Muller, 

1764) and not to C. California Westwood, [1851] (Fig. 14). These populations are currently attributed to 
subspecies C. tullia kodiakVJ . H. Edwards, 1869 (type locality USA: Alaska, Kodiak), C. tullia mixturata 
Alpheraky, 1897 (type locality Russia: Kamchatka; it remains to be investigated if this name applies to 
the Nearctic taxon) and C. tullia yukonensis W. Holland, 1900 (type locality Canada: Yukon and USA: 
Alaska). Compared to European C. tullia, they show Fst/Gmin of 0.14/0.05, which are in the range for 
conspccific populations, however with a more limited gene exchange than typical (Gmin is less than 0.1). 
Therefore, until further research shows otherwise, we leave these three subspecies with C. tullia. Future 
genomic studies of Coenonympha tullia viluiensis Menetries, 1859 (type locality Russia: Vilyuy River) 
are also needed, because if it falls in the same clade with the three subspecies of C. tullia from North 
America, it would be the oldest name in this clade. And if this clade is found to be a species distinct from 
C. tullia, it will be the name of this species. In summary, according to our genomic analysis, both C. tullia 
and C. California are present in North America, and all American taxa of the tullia complex other than the 
northernmost subspecies C. t. kodiak, C. t. mixturata and C. t. yukonensis belong to C. California. 


i— | 0.9 

Ag= 


tympha tullia tullia119097E08|M[Sweden,Lockhyttan|1978 
jenonympha tullia tullial 1909/E09|M|France| 1975 
Coenonympha tullia tullia] 19097E10|M|Czech Republic,Bohemial 1976 

Coenonympha tullia kodiakl 19098H03|USA:AK,Kodiak Island] 1899 





- Coenonympha tullia mixturata?|i9098H05|USA:AK,Kotzebue|195i 

- Coenonympha tullia yukonensis] 19098H06|USA:Canada,Yukon Terr.11961 
Coenonympha tullia yukonensis|19098D12|Canada:Yukon Terr.|1961 

Coenonympha tullia California119097E12jF|USA:CA,Marin Co.11993 
Coenonympha tullia ampelosj 19098HlljUSA:CA.PIumas Co.11984 
Coenonympha tullia yontocket|19098H09|PT|USA:CA,Del Norte Co.| 1979|USNM 
Coenonympha tullia eunomia| 19098H07jUSA:OR,Lane Co.11994 
Coenonympha tullia ochracea|19098Dli|USA:MT,Ravalli Co.11981 
Coenonympha tullia ochracea[9388|F|USA:WY,Park Co.,Yellowstone National Park|2017 
Coenonympha tullia beniamini|19098G12|M[Canada:AlbertaJ2003 
‘ L " " ilfia t — : — - - - - J - - 1 ” 


Coenonympha tullia benjamini|19098bl0|M|Canada:Saskatchewan| 1958 


Coenonympha tullia benjaminiJ19098H01JM|Canada:Alberta|1991 
Coenonympha tullia bemamini[19098D08[USA:NE,Cherry Co.11985 


iiaa- 


Coenonympha tullia inornata 19098D05 |USa:MI,P resque Isle Co.J 1978 
Coenonympha tullia inornata] 19098D03|USA:NH.Mt. Washington! 1982 

Coenonympha.tullia inornatal19098D02|USA:VT,Addison Co.1 1974 


Coenonympha tullia inornatal 19098D04|USA:NY,Niagra Co.|2005 
- Coenonympha tullia mcisaacil 19098C07|Canada:Newfoundland| 1975 

- Coenonympha tullia mcisaaci[19098C09[Canada:Newfoundland|1971 

- Coenonympha tullia nipisiquit]01506C07|HT|M|Canada:New Brunswick,Bathurst| 1939 

- Coenonympha tullia nipisiquit[19098D01|M|Canada:New Brunswickl 1952 

Coenonympha tullia nipisiquitl 19098C12|M|Canada:New Brunswickl 1954 

Fig. 14. Coenonympha tullia (blue and magenta) and California (red). 
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Palaeonympha Butler, 1871 is a subgenus of Megisto Hiibner, [1819] 

The New World genus Megisto Hiibner, [1819] (type species Papilio eutytus Fabricius, 1775, a junior 
homonym, considered a synonym of Papilio cvmela Cramer, 1777) became monotypic after the transfer 
of Euptychia ruhricata W. H. Edwards, 1871 
(type locality USA: Texas, McFennan Co.) to 
Cissia Doubleday, 1848 (type species Papilio 
clarissa Cramer, [1780], a junior subjective 
synonym of Papilio penelope Fabricius, 1775) 

(Zacca et al. 2018). A genus Palaeonympha 
Butler, 1871 was proposed for a newly described 
species P. opalina Butler, 1871 from China and 
remained monotypic since. In wing patterns, it 



Palaeonympha opalina 119114B02[Taiwan11957 
Palaeonympha opalina|19114B03|Taiwan|1957 

Megisto cvmela cvmela 43501 USA: IN .Lake Co.12015 
[Megisto cymela cymelal853lluSA:TX,Harrison Co.|2017 
Megisto cymela viola[8622|USA:FL,Levy Co.|2017 

- Vanima labe[SAMN0974-5406]Ecuador:Carchi,Lita|2014 
Vanima palladia|SAMN0974544-3[Brazil:Sao Paulo|2013 
Yphthimoides pacta|SAMN09745,451 [Brazil:Mato Grosso12013 
Yphthimoides renata|19117G08 Costa Rica:ACG|2011 
Yphthimoides borasta|SAMN0974-5452|Brazil:Minas Gerais|2014 
Cissia cleophes[19118B01]Me>cico:Guerrero|1906 
Cissia proba|SAMN09745391]Ecuadorj2014 
Cissia rubricata rubricataj8382 USA:TX,Blanco Co.|2017 

- Cissia penelope|SAMN09745395|Ecuador|2014 

- C i s s i a e o u s | BAM 11 'j 9 ~45454| B ra ill: S a o P a u I o 12 014 

Hermeuptychia sosybius[3310|USA:TXJ>abine Co.|2015 
Hermeuptychia intricataliseslPTlUSAifx^ort Bend Co.|2013 

Fig. 15. Palaeonympha (red) is Megisto (blue); Cissia (magen¬ 
ta), Vanima (brown). IDs with "SA" are from NCBI database. 


22 















































bears an uncanny resemblance to M. cymela (Fig. 16), and was indeed sister to Megisto both in gene 
marker-based (Zacca et al. 2018) and genome-scale phylogenies (Espeland et al. 2019). Genomic 
sequencing confirms a close relationship between monotypic Megisto from the New World and monotypic 
Palaeonvmpha from the Old World (Fig. 15 blue and red), closer than many species of Cissia are to each 
other (Zacca et al. 2018; 

Espeland et al. 2019) and 
about the same divergence 
as two species of Vanima 
Zacca, Casagrande & O. 

Mielke (Fig. 15 brown). 

The COI barcodes of 
cymela and opalina differ 
by 8.8% (58 bp), within 
the range of many con¬ 
geners. To emphasize this 
close kinship between 
these two disjunct species 
(i cymela and opalina) 
apparent from their geno¬ 
types and phenotypes, we 
propose that, instead of 
each being placed in its 
own monotypic genus, they are congeneric. We assign a new status of subgenus to Palaeonympha, 
resulting in Megisto opalina Butler, 1871, new combination. Thus, Megisto becomes a Holarctic genus, 
yet again indicating elaborate connections between the Old and the New Worlds. 



Megisto (Megisto) cymela M. (P.) opalina 


Fig. 16. Megisto opalina (Taiwan, sequenced as NVG-19114B02 & NVG-19114B03) 
_ and M. cymela (USA: OK, Carter Co., Murray Lake, 9-Apr-1999). _ 


Cissia cleophes (Godman & Salvin, 1889) does not belong to Megisto Hiibner, [1819] 



Placed in Megisto recently on the basis of morphological similarities (Zacca et al. 2020), Euptychia 
cleophes Godman & Salvin, 1889 (type locality Mexico, Guerrero, Fig. 17) is not monophyletic with the 
type species of Megisto (Fig. 15 
magenta vs. blue) and therefore 
does not belong to Megisto, 
unless this genus is expanded to 
include species placed in other 
genera, such as Yphthimoides 
Forster, 1964, Vanima and 
Cissia. In our tree (Fig. 15), 
cleophes falls within Cissia : it is 
in the clade with C. penelope, 
which is a valid name for the type 
species of Cissia. Therefore, because the tree strongly supports that cleophes does not belong to Megisto, 
we suggest to leave this species in Cissia (where it may belong and was placed previously, but statistical 
support for this clade is weaker) pending genomic analysis of other species from this and nearby clades. 


Fig. 17. Cissia cleophes sequenced as NVG-19118B01. 


Hyponephele Muschamp, 1915 is a subgenus of Cercyonis Scudder, 1875 

In the genomic tree (Fig. 18), Hyponephele Muschamp, 1915 (type species Papilio lycaon Rottemburg, 
1775) is sister to Cercyonis Scudder, 1875 (type species Papilio alope Fabricius, 1793, placed as a 
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subspecies of Papilio pegala Fabricius, 1775), in agreement with previous findings (Pena et al. 2006). 
The tree reveals that the genetic divergence between Hyponephele and Cercyonis is smaller than that 
within Erebia Dalman, 1816 (type species Papilio 
ligea Linnaeus, 1758). The COI barcodes of//, lycaon 
and C. pegala differ by 7.9% (52 bp). We give 
Hyponephele a new status of a subgenus within 
Cercyonis. This change may not be welcomed by the 
Old World Lepidopterists who are used to the name 
Hyponephele applied to its many species, similar to 
how Argynnis is not welcomed in America to include 
Speyeria as its subgenus. However, this name change 
highlights the close relationship between the two 
subgenera ( Hyponephele and Cercyonis) making 
Cercyonis a Holarctic genus, similar to Erebia 
Dalman, 1816 (type species Papilio ligea Linnaeus, 1758) in divergence and distribution. This is yet 
another step towards more internally consistent genus-level classification in butterflies. 


Cercyonis pegala texana|4506|USA:TX,Wise Co.|2015 
Cercyonis oetus oetus|PA0112|USA:NV,Elko Co.|2016 
"ercyonis incognita|PA0435|USA:CA,Mendocino Co.|2017 
rcyonis sthenele|7259|ST|USA:CA / San Fransico Co.|old 
Jercyonis meadii|9708|USA:TX Jeff Davis Co.|2017 
Hyponephele lycaonl 19113H06|North Macedonia|1980 
Hyponephele dysdoral 19115E11 (Tajikistan11975 
Maniola telmessial 19115F05|Greece11983 
Maniola jurtina|20035B06|France:Var|2018 
Pyronia tithonus| 19113H10|France|2008 
Pyronia bathseba|19113H09|Spain|1976 

Pyronia cecilia 119113H08| Spain 11976 

Erebia callias|9376|USA:WY,Park Co.|2017 
Erebia vidleri|8001|USA:WA,Whatcom Co.11980 
Erebia occulta|16106G09]Canada:Yukon Territory12016 
• Erebia epipsodea brucei|6457|USA:CO,Grand Co.|2016 
Erebia magdalena|6549|USA:CO / Park Co.|2016 
Erebia discoidalis|7999|USA:WI,Oneida Co.11988 
Erebia rossiil 16106E05|Canada:Yukon Territory12016 
Erebia disal 16106E08|Canada:Yukon Territory]2016 

Fig. 18. Cercyonis (blue), Hyponephele (red), Maniola 
(green), Pyronia (magenta), and Aphantopus (cyan). 



Pyronia Hiibner, [1819] and Aphantopus Wallengren, 1853 
are subgenera of Maniola Schrank, 1801 

The genomic tree reveals that Pyronia Hiibner, [1819] (type species Papilio tithonus Linnaeus, 1771) is 
not monophyletic (Fig. 18), and Pyronia bathseba (Fabricius, 1793) is sister to Aphantopus Wallengren, 
1853 (type species Papilio hyperantus Linnaeus, 1758) with strong support (Fig. 18). Two other Pyronia 
species we sequenced are not monophyletic either: the type species of the genus is sister to Maniola 
Schrank, 1801 (type species Maniola lemur Schrank, 1801, which is a junior subjective synonym of 
Papilio jurtina Linnaeus, 1758). To restore monophyly, it is possible to break Pyronia into smaller 
genera, and these already have names available: Pasiphana de Lesse, 1952 (type species Papilio bathseba 
Fabricius, 1793) and Idata de Lesse, 1952 (type species Epinephele ida var. cecilia Vallantin, 1894). 
Alternatively, they can be grouped in some ways to form more inclusive monophyletic genera. The tree 
(Fig. 18) reveals three clusters of species of equivalent rank. One of these clusters is the genus Erebia. 
The other one is the genus Cercyonis as we presently define it (including Hyponephele as a subgenus). 
Therefore, it is meaningful to treat the third group as a single genus as well. Hence, we propose that 
Aphantopus, Pyronia, Pasiphana and Idata are subgenera of Maniola, new status. It is unfortunate that 
genomic data suggest abandoning the familiar Aphantopus and Pyronia as genera, however breaking this 
more inclusive but genetically prominent genus Maniola into four or five very small genera is even less 
appealing to us. 


Dione incarnata N. Riley, 1926 is a species 
distinct from Dione vanillae (Linnaeus, 1758) 


Inspection of a nuclear genomic tree reveals pronounced divergence between Dione vanillae (Linnaeus, 
1758) (type locality South 
America, probably Surinam) 
specimens from Jamaica and 
from the USA. (Fig. 19). This 
divergence is of about the 
same magnitude as that 
between the distinct species 
Dione moneta Hiibner, [1825] 
and Dione juno (Cramer, 1779), and somewhat smaller than that between more distant relatives 
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I r 1 .Heliconius charithonia tuckeril 46951 USA: FL,Levy Co,Yankeetown|2015 

IpT) -Heliconius charithonia tuckeri|8186|M|USA:FL,Miami-Dade Co.,4.8 air mi E of Florida|2017 

— 1 Ir-o ^liconius charithonia vazquezae|4474|M|USA:TX,Travis Co,Austin|2015 

‘—'Heliconius charithonia vazquezae|5097|M|USA:TX,Starr Co,Rio Grande City|2015 
- Heliconius erato petiverana|17117B05|M|Mexico:San Luis Potosi|1976 

— - Dione moneta} 14114D10|USA:TX,Hidalgo Co.,Mission|2014 

- Dione juno huascumal 17117B01|M|Mexico:San Luis Potosi|1978 

|—, Dione vanillae insularis|10475|Jamaica:12|2017 
I Ho §i one van ‘H ae insularis|10272|Jamaica:2|2017 

. *— 'Dione vanillae insularisj 10536|Jamaica:3|2017 

I j o <Q[one vanillae nigrior[8139lF[USA:FL,Liberty Co.,SH12|2017 

|_ y-r Dione vanillae nigrior|8175|F|USA:FL,Miami-Dade Co.,6 mi SW of Florida City12017 

[i—= -Dione vanillae incarnata| 10145|M|USA:AZ,Santa Cruz Co.,Pajarito Mtns.|2017 
1 —•' Dione vanillae incarnata1345l|M|USA:TX,Hidalgo Co,Old Rio Rico Rd. 1.5 air mi|2015 

Fig. 19. Dione incarnata (red) and D. vanillae insularis (blue). 




























Heliconius charitonia (Linnaeus, 1767) and //. erato (Linnaeus, 1758) (Fig. 19). In contrast to their 
separation from Dione vanillae insularis (Maynard, 1889) (type locality Bahamas, Fig. 19 blue), the name 
currently applied to Jamaican populations, the two USA subspecies: eastern Dione vanillae nigrior 
(Michener, 1942) (type locality USA: Florida, Monroe Co.) and western Dione vanillae incarnata N. 
Riley, 1926 (type locality Mexico: Durango), cluster closely with each other (Fig. 19 red). The Fst/Gmin 
statistics for the two groups (red vs. blue) are 0.81/0.0002, indicating nearly absent gene exchange 
between these groups, the smallest of all sister species we studied in this work. Search of the BOLD 
database (Ratnasingham and Hebert 2007) reveals that Jamaican specimens group closely by their COI 
barcodes with specimens from other Caribbean Islands (Cuba, Dominican Republic, Puerto Rico), and in 
particular with specimens from the Bahamas (the type locality of D. v. insularis. Therefore, Jamaican 
populations may indeed be referred to as D. v. insularis. The COI barcode differences between Jamaican 
insularis and USA incarnata is 2.9% (19 bp). For these reasons, D. vanillae insularis is apparently a 
species distinct from the USA species consisting of subspecies nigrior and incarnata. 

These two species ( insularis and the USA species) are distinct phenotypically with wing shapes and 
patterns, as described by Maynard (1889). In fact, Maynard proposed Agraulis insularis as a species, and 
he considered the USA species to be Agraulis vanillae, characterized by longer wings, smaller black 
spots, white dots in forewing black spots compared to insularis. However, both the Merian illustration 
(1705), which was mentioned in the original description of D. vanillae (Linnaeus, 1758), and the 
lectotype specimen (Honey and Scoble 2001) reveal the insularis phenotype (e.g. extended black band on 
the forewing, broader wings) that differs from the USA specimens. In his revision, Michener gave a key 
to the vanillae complex taxa (Michener 1942). The first doublet separates incarnata with nigrior from all 
other taxa, including the nominotypical vanillae. Therefore, the USA specimens are not D. vanillae, but a 
species that may be referred to by the oldest name applicable to North American populations: D. 
incarnata, and we currently leave insularis as a subspecies of D. vanillae. Michener also raised the 
possibility that incarnata may be a distinct species, based on incarnata specimen from Colombia, in 
which case it would be sympatric with D. vanillae. However, these Colombian specimen might have been 
mislabeled. Due to the evidence presented here, Dione incarnata N. Riley, 1926, new status, appears to 
be a distinct species, with D. incarnata nigrior (Michener, 1942), new combination, being its subspecies. 


Danaus eresimus (Cramer, 1777) belongs to subgenus Danaus Kluk, 1780, together 
with Danausplexippus (Linnaeus, 1758), and not to subgenus Anosia Hiibner, 1816 

together with Danaus gilippus (Cramer, 1775) 

Danaus eresimus (Cramer, 1777) (type locality Suriname) looks superficially similar to Danaus gilippus 
(Cramer, 1775) (type locality Brazil: Rio de Janeiro) and at times it is a challenge to distinguish these two 
species. On the contrary, Danaus 
plexippus (Linnaeus, 1758) (type 
locality USA: New York, Orleans Co.) 
is superficially more different from 
either D. eresimus or D. gilippus 
(Warren et al. 2016). Due to these 
superficial similarities and differences, 
traditionally, only the former, as the type species, belonged to Danaus Kluk, 1780, and the latter two 
species were placed in Anosia Hiibner, 1816 (type species Papilio gilippus Cramer, 1775) (Ackery and 
Vane-Wright 1984; Pelham 2008; Pelham 2020). However, among these three species, genomic data 
(both nuclear and mitochondrial genomes) place D. eresimus as a sister to D. plexippus with high 
confidence (Fig. 20), and D. gilippus is a sister to that clade of the two species, in agreement with 
previous DNA-based analyses (Zhan et al. 2014; Aardema and Andolfatto 2016). Therefore, we transfer 
Danaus eresimus (Cramer, 1777) from the subgenus Anosia where it does not belong, to the subgenus 
Danaus in accord with the phylogeny of these three species. This change has already been implemented 
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by Pelham in the most recent version of the catalogue (2020) after the discussion of our genomic data and 
previous works with NVG. Here, we simply formalize this change in a publication. 


Chlosyne coronado (M. Smith & Brock, 1988) is a species 
distinct from Chlosyne fulvia (W. H. Edwards, 1879) 


fiulvia and coronado are 
0.49/0.015, comparable 
to those of fulvia vs. 
leanira : 0.40/0.012. COI 

barcodes are 2.7% (18 bp) different between fulvia and coronado, but are 0% between coronado and 
pariaensis. Furthermore, according to the map of sequenced specimens (Fig. 21 ), fulvia and coronado 
may be sympatric in north-central Arizona. For all of these reasons, we propose species status for 
Chlosyne coronado (M. Smith & Brock, 1988) new status, and Chlosyne coronado pariaensis (M. Smith 
& Brock, 1988) new combination, as its subspecies. 




AZ: Coconino Co. 


NVG-20011F12 


AZ: Mohave Co., Peach Springs 


leanira 


pariaensis 


almai 


20011F12 


-pelona 


20011G01 


i20011Fll. 


wrightii 


20011G02' 


coronadoi 


Chlosyne coronado 


Chlosyne fulvia 


Chlosyne 

taxa 

| - wrightii 
□ - elegans 
2 - leanira 
2 - alma 
/\ -pariaensis 
f\ - coronado 
O - fulvia 


Chlosyne leanira wrightii118034C09|CA:San Diego Co. 
Chlosyne leanira wrightii (=pelonaJ 19061Bll|HTlCA:LACo. 
-u,--' 19061H09[AT|CA:SLO Co. 


Chlosyne leanira € 


In the leanira group, Chlosyne fulvia coronado (M. Smith & Brock, 1988) (type locality USA: Arizona, 
Pima Co.) has been a 
subspecies-level taxon 
since its description 
(Smith and Brock 1988). 

The genomic tree of the 
group revealed that 
Chlosyne fulvia (W. H. 

Edwards, 1879) (type 
locality USA: “Western 
Texas”) is paraphyletic 
with respect to Chlosyne 
cyneas (Godman & 

Salvin, 1878) (type 
locality Mexico: Oaxaca), 
and C. fulvia coronado 
with Chlosyne fulvia 
pariaensis (M. Smith & 

Brock, 1988) (type 
locality USA: Utah, 

Kane Co.) form a clade 
more prominent than 
Chlosyne leanira (C. 

Felder & R. Felder, 

1860) (type locality 
USA: California, Plumas 
Co.) (Fig. 21). The Fst/ 

Gmin statistics for 


, elegans]. 

Chlosyne leanira leanira|PA0352JCA:Stanislaus Co. 
Chlosyne leanira alma| 18034D04|NV:Esmeralda Co. 
Chlosyne leanira alma| 18034D05]NV:Nye Co. 

Chlosyne fulvia pariaensis| 17109D03|HTJUT:Kane Co. 

Chlosyne fulvia coronado|20011F12|AZ:Mohave Co. 

— Chlosyne fulvia coronado] 20011 FI 11 AZ: Coconino Co. 
Chlosyne fulvia coronado|_17109D02|HT|AZ:Pima Co. 

— Chlosyne fulvia fulvia|20011G01|AZ:Coconino Co. 
Chlosyne fulvia fulviaI20011G02IAZ:Apache Co. 
Chlosyne fulvia fulvia117117D03|TX:Hall Co. 

Chlosyne fulvia fulvia 117117D02JTX:Brewster Co. 
Chlosyne fulvia fulvia 110956|TX: Brewster Co. 
Chlosyne fulvia fulvia 112371jKS: Barber Co. 

Chlosyne fulvia fulvia117117D04ITX:Terrell Co. 

- Chlosyne fulvia fulviafl7117D05|TX:Terrell Co. 
Chlosyne cyneas|17117D011Mexico:Tlaxcala 


fulvia. 

Western Texas 


o - type locality: 

specimens sequenced 
o - type locality: 

no specimens sequenced 

type localities are indicated 
by tiny circles placed inside 
specimen symbols, or small 
circles framed with color of taxa 
(if no specimens from these 
localities were sequenced), labeled 
with names of taxa preceded by = 
for junior subjective synonyms, 
shown in smaller, thinner font 

Fig. 21. Chlosyne fulvia and coronado with relatives. Four specimens from Arizona are 
shown (dorsal: right and ventral: left) and are labeled with their DNA sample numbers on the 
map that shows the localities of all sequenced specimens from the USA included in the tree. 


Chlosyne chinatiensis (Tinkham, 1944) is a species 
distinct from Chlosyne theona (Menetries, 1855) 

In their revision of the theona group, Austin and Smith (Austin and Smith 1998) placed Melitaea 
chinatiensis Tinkham, 1944 (type locality USA: Texas, Presidio Co.) as a subspecies of Chlosyne theona 
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(Menetries, 1855) (type locality Nicaragua). While accepting subspecies-level treatment of chinatiensis, 
Pelham (2008) writes: "There is considerable 
reason to consider this a distinct species from 
theona. More investigation is required." We 
carried out our genomic investigation by 
sequencing of three chinatiensis specimens from 
the US and Mexico (Fig. 22) and found that their comparison with theona specimens from across the 
range (from Arizona, Texas and Costa Rica) results in the following Fst/Gmin statistics: 0.35/0.019, 
indicating genetic differentiation and low gene exchange consistent with chinatiensis being a species-level 
taxon. The COI barcodes of C. chinatiensis and C. theona thekla differ by 1.4% (9 bp), but those of C. 
chinatiensis and C. theona bolli differ by 0.6% (4 bp). Moreover, C. chinatiensis is sympatric with C. 
theona bolli in west Texas, e.g. in the Big Bend National Park. Given this evidence, we reinstate 
Chlosyne chinatiensis (Tinkham, 1944) as a species. 



Family Hesperiidae Latreille, 1809 

Phocides lilea (Reakirt, [1867]) is a species 
distinct from Phocides polybius (Fabricius, 1793) 

A tree constructed from protein-coding regions of their nuclear genomes reveals that specimens of 
Phocides polybius (Fabricius, 1793) (type locality "Indiis", likely Suriname) partition into two distinct 
clades (Fig. 23). One of the clades 
consists of Phocides polybius lilea 
(Reakirt, [1867]) (type locality 
Mexico: Veracruz, Fig. 23 red). The 
other clade includes all other taxa of 
this species (Fig. 23 blue). Fst/Gmin 
statistics for the Z chromosome 
comparison of these two clades are 
0.39/0.021, suggesting that they 
represent two distinct species. 

Moreover, COI barcodes of lilea from Mexico and polybius from Guyana show 3% (20 bp) difference. 
Therefore, we conclude that Phocides lilea (Reakirt, [1867]) is a species-level taxon, reinstated status. 

Furthermore, we sequenced the only known syntype of an enigmatic taxon Eiycides imbreus Plotz, 
1879 from the ZMFIB collection, illustrated in Warren et al. (2016). This specimen is a syntype because it 
is a uniquely patterned specimen that carries appropriate labels, agrees with the original description and 
looks similar to the unpublished Godman copy of the Plotz illustration (in BMNH, inspected by NVG) 
(Godman 1907). It is an unusual specimen lacking an orange bar in the forewing discal cell (usually 
extending to costa) typical for P. polybius. Evans (1952) treated this name as a distinct species Phocides 
imbreus Plotz, 1879 based on a rather poor illustration of this specimen in Draudt (1921)—it is unlikely 
that Evans saw the actual specimen. Mielke & Casagrande (2002) inspected the syntype and synonymized 
the name with P. polybius lilea due to general phenotypic similarity and the lack of orange coloring on the 
fringe around the hindwing tornus. According to our genomic results (Fig. 23), imbreus is confidently 
placed with Phocides polybius polybius, revised placement of a synonym, and is probably an aberrant 
specimen of polybius, not lilea, lacking any orange coloration on its wings, not just on the fringe, but also 
a forewing orange bar. However, the head of the syntype retains the usual orange patterns including 
orange palpi and cheeks. Our revised synonymy is further supported by the label data on the specimen 
stating "Am. m.", which probably stands for America meridionalis (South America), where P. lilea is not 
known to occur. 



Phocides polybius lilea (=sanguinea)|15096C08|T|USA:TXlprior to 1872JCMNH 
— Phocides polybius lilea r=albiciliata)J18094C09|HT|F[Guatemala|oldJMTD 
Phocides polybius lilea (=albicilla)115029B121ST|Mexico]prior to 1869|ZMHB 
- Phocides polybius lilea| 17113G02lF)USA:TX,Cameron Co.11986 

Phocides polybius polybius|19087H10|M|Venezuela,Merida|old 
— Phocides polybius polybius|17098G01|Panama| 1978 
Phocides polybius polybius[19087H08|MJFrench Guiana|old 
— Phocides polybius polybius|19087Hli|M|French Guiana,Cayenne|old 
- Phocides polybius polybius (x=gunderi)| 15097B03|T|Guyana|1927|CMNH 
Phocides polybius polybius11908/H09|Guyana|old 
Phocides polybius polybius (=imbreus)|15029B09|ST|"South America"|ZMHB 
Phocides polybius phanias|17098G02|Brazil:Sao Paulo|1991 
■fU Phocides polybius phanias (=unimacula)| 18086F05|?ST|M|Argentina| 1906IMNHP 
— Phocides polybius phanias (=unimacula)| 18086F07|?ST|M|Argentina| 1906|MNHP 
Phocides distans licinus117099D07|Panama| 1969 
Phocides pialia pialia (=zancleius)| 18025C09|HT|M|Brazil:SC|before 1932|AMNH 
Phocides yokhara|18031C04|Ecuaaor|2002 
- Phocides perillus|18039D06|M|Colombia,Bogota|old 

Fig. 23. Phocides lilea (red) and polybius (blue). 
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Cecropterus nevada (Scudder, 1872) and Cecropterus dobra (Evans, 1952) 
are species distinct from Cecropterus mexicana (Herrich-Schaffer, 1869) 


Genomic analysis of Cecropterus mexicana (Herrich-Schaffer, 1869) (type locality Mexico) reveals a 
pronounced divergence between its subspecies (Fig. 24) that was analyzed further. The genomic tree 
shows separation between some 
of them comparable to that from 
Cecropterus diversus (E. Bell, 

1927) (type locality USA: 

California, Plumas Co.). While 
Cecropterus mexicana aemilea 
(Skinner, 1893) (type locality 
USA: Oregon, Klamath Co., 
male syntype sequenced), C. mexicana blanca (J. Scott, 1981) (type locality USA: California, Mono Co.) 
and C. mexicana nevada (Scudder, 1872) (type locality USA: California, Sierra Nevada) group closely 
together (all three unified under the name nevada below), C. mexicana dobra (Evans, 1952) (type locality 
USA: Arizona, Graham Co.) forms a clade distinct from them and C. mexicana. The Fst/Gmin statistics 
for these clades are: mexicana vs. dobra : 0.34/0.021, mexicana vs. nevada : 0.37/0.010, nevada vs. dobra : 
0.30/0.055. We see that nevada and dobra exchange genes more frequently with each other than do each 
of them with mexicana. Differences between COI barcodes in pairs of these species are: mexicana and 
dobra : 1.8% (12 bp), mexicana and nevada : 1.1% (7 bp), nevada and dobra : 1.7% (11 bp). For 
comparison, the COI barcodes of aemilea, blanca and nevada are 100% identical. Curiously, in contrast 
to nuclear genomes (Fig. 24), mitochondrial genomes (as reflected by barcodes) place mexicana closer to 
nevada, and dobra farther away from them, which is yet another example of the peculiarity of 
mitochondrial evolution. Deriving further support from genitalic and wing pattern differences mentioned 
by Evans (1952), we suggest that Cecropterus nevada (Scudder, 1872), reinstated status, and 
Cecropterus dobra (Evans, 1952), new status, are species-level taxa, not subspecies of Cecropterus 
mexicana (Herrich-Schaffer, 1869). Then, we treat Cecropterus nevada aemilea (Skinner, 1893) and 
Cecropterus nevada blanca (J. Scott, 1981), new combinations, as subspecies of C. nevada. 


Cecropterus mexicana aemilea|15095F02|ST|USA:OR,Klamath Co.|before 1893|CMNH 
Cecropterus mexicana blanca|PA01151|USA:CA,Mono Co..White Mountains|2019 
Cecropterus mexicana blanca|PA01152|USA:CA,Mono Co..White Mountains|2019 
Cecropterus mexicana nevada|14114A09 FJUSAiCA,Tulare Co..Sherman Pass area|2003 
Cecropterus mexicana nevada 14114A10 M|USA:CA,Tulare Co..Sherman Pass area|2003 
Cecropterus mexicana dobra114114A07|M USA:AZ,Coconino Co.,Lockett Label 1993 
- Cecropterus mexicana dobra|11182|M|USA:NM,Otero Co..Lincoln National Forest|2018 

- Cecropterus mexicana mexicana (=ananius]| 15032A01|ST|M|Mexico|oldlZMHB 

Cecropterus mexicana mexicana|14114B01IM|Mexico:Oaxaca,San lose Pacifico] 1986 
Cecropterus mexicana mexicana|14114A12|Mexico;Flidalgo,Cuesta Colorada|1981 
Cecropterus mexicana mexicana|14114All)F|Mexico:F)idalgoJacala 

'-...USAiCA,Madera Co.|2007 ~ 

- “ .McGee Creek|2015 

ulla Co.,3 mi N of Panaceal 1980 

Fig. 24. Cecropterus nevada (red), dobra (magenta) and mexicana (blue). 



Telegonus anausis Godman & Salvin, 1896, is a species 
distinct from Telegonus anaphus (Cramer, 1777), and 
Telegonus anausis annetta (Evans, 1952) is its subspecies 


A polytypic species Telegonus anaphus 
includes four subspecies (Evans 1952) 
single representative of the USA 
fauna: Astraptes anaphus annetta 
Evans, 1952 (type locality Costa 
Rica). The genomic tree of these 
taxa reveals a split between the 
nominotypical T. anaphus and all 
others (Fig. 25 blue vs. red clade). 
The oldest name in the second (red) 
clade is Telegonus anausis Godman 
& Salvin, 1896 (type locality St. 
Vincent, Grenada, Dominica, His- 


(Cramer, 1777) (type locality Suriname), in addition to nominal, 
all of which we sequenced and analyzed. Among them, there is a 


°80.14 

D.W 


Telegonus anaphus annetta] 14111H04)M|USA:TX,Starr Co.,Ft. Ringgold|2014 
Telegonus anaphus annetta|14111H08(Mexico:Tamaulipas 
Telegonus anaphus annetta|19075D10|Honduras|1979 


Telegonus anaphus annetta|18121A02|Costa Rica:ACG|2007|07-SRNP-20089 
Telegonus anaphus annettal 17106C09|Costa Rica:ACG|2010| 10-SRNP-67024 
Telegonus anaphus annetta|19075D12|Colombia|1992 
Telegonus anaphus annettal 14113F08|Mexico:Guerrero| 1985 




— Telegonus anaphus anoma| 19075E03|Trinidad| 1982 

-- Telegonus anaphus anoma|MJWC15-52|M|Trinidad,Point Gourdel 1993 

I— 5 - 5 ^ Telegonus anaphus anausis|19075F06|Puerto Rico|1934 

—P-H- Telegonus anaphus anausis[19075F04|Cuba,Santiago|old 

y4s- Telegonus anaphus anausis|10688|F|Jamaica,Mandeville|2018 

Telegonus anaphus anausis|19075F07|Grenada,Carriacou|1921 


- Telegonus anaphus aniza|18028H01|Ecuador|1988 

-- Telegonus anaphus aniza|19075E09|Peru|1999 

TTi- Telegonus anaphus anaphus| 19075Ell|Brazil:Paraiba| 1956 

X72- Telegonus anaphus (=anaphides)|18086G03|HT|BR:Bahia|1912|MNHP 

n- Telegonus anaphus anaphus|19075E10|Brazil:Rondonia|1993 

1 - Telegonus anaphus anaphusl 19087A01|French Guiana.,Cayenne|old 

- Telegonus jaira|10323|M|Jamaica|2017 

Fig. 25. Telegonus anausis subspecies (red branches) and anaphus (blue). 


paniola). The Fst/Gmin statistics computed on some of these subspecies are: anaphus vs. annetta 
0.40/0.01 and anaphus vs. anausis 0.45/0.009, indicating species-level differentiation. The COI barcodes 
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are 4.3% (28 bp) different between anaphus and annetta, which is larger than a typical difference between 
closely related species. Comparing subspecies within the second cluster (Fig. 25 red clade), we were not 
able to gather evidence for their change of status to species. For instance, Fst/Gmin of anausis vs. annetta 
is 0.16/0.11, which by itself does not justify Central American annetta as a species distinct from the 
Caribbean anausis. Therefore, we reinstate Telegonus anausis Godman & Salvin, 1896 as a species and 
transfer all (except the nominotypical) subspecies currently placed in T. anaphus to T. anausis to fonn the 
following new combinations: Telegonus anausis annetta (Evans, 1952), Telegonus anausis anoma 
(Evans, 1952), and Telegonus anausis aniza (Evans, 1952). A curious accident here is that as a result of 
genomic work, the name originally proposed (Godman and Salvin 1896) as Telegonus anausis is now 
returned to its original combination and status nearly 125 years later. 


Epargyreus huachuca Dixon, 1955 is a species 
distinct from Epargyreus clarus (Cramer, 1775) 

We obtained whole genome shotgun sequences of specimens from all known distinct groups of US 
populations currently assigned to Epargyreus clarus (Cramer, 1775) (type locality "Suriname", later 
corrected to USA: Virginia, 

Rockingham Co.). A tree made 
from protein-coding regions of 
the Z chromosome (Fig. 26) 
reveals a prominent split into 
two clades (blue and red) with 
divergence comparable to that 
between Epargyreus orizaba 
Scudder, 1872 and Epargyreus 
cruza Evans, 1952, shown as 
outgroups. The holotypes of 
both E. clarus huachuca Dixon, 

1955 (type locality USA: Arizona, Cochise Co.) and E. clarus profugus Austin, 1998 (type locality USA: 
Nevada, Clark Co.) have been sequenced, along with a possible type specimen of E. clarus californicus 
MacNeill, 1975 (type locality USA: California, El Dorado Co.). The Fst and Gmin between E. clarus 
clarus and E. clarus huachuca are, 0.35 and 0.03 respectively, compared to those of 0.04 and 0.16 
between E. clarus huachuca and E. clarus profugus. The difference in COI barcodes of E. clarus clarus 
and E. clarus huachuca is about 2% (~ 13 bp difference). For these reasons, we suggest that Epargyreus 
huachuca Dixon, 1955 is a distinct species, new status, and it includes Epargyreus huachuca profugus 
Austin, 1998, new combination, as a subspecies. We found that genetic differentiation between 
populations of E. clarus californicus is lower than between others (Fig. 26, a tight cluster of specimens 
with shorter terminal branches) suggesting either a recent dispersal or a bottleneck. Whether this genetic 
purge resulted in any degree of reproductive isolation of E. clarus californicus from other groups of E. 
clarus populations remains to be investigated. 


Epargyreus clarus califomicusl 18037E12|USA:CA,Riverside Co. 


u r Epargyreus clarus califomicusl 15117A04|USA:CA,Tulare Co. 

- Epargyreus clarus califomicusl 15105E06|T?|USA:CA,Eldorado Co 



Epargyreus clarus califomicusl 15116Hll|USA:CA,Tulare Co. 

- Epargyreus clarus|19125A01|Canada:Ontario 

Epargyreus clarus clarus|4283|USA:IN,Montgomery Co. 

Epargyreus clarus clarus|4192|USA:Tx,Dallas Co, 

Epargyreus clarus clarus|3526|USA:TX,San Jacinto Co. 

Epargyreus clarus clarus161551USA:VA ; Rockbridge Co. 
Epargyreus clarus clarus|9230|USA:AR,Polk Co. 

Epargyreus clarusl 19124H12 USA:MT ; Carbon Co. 
Epargyreus clarusf 19039H02[USA:MT,Lincoln Co. 

Epargyreus clarus| 19039H03|USA:CO,Gunnison Co. 
Epargyreus clarusl 19039H04|USA:CO.Boulder Co. 

Epargyreus clarus huachuca|_15095E07|HT|USA:AZ t Cochise Co. 
- . -,- L - , 15105E08|USA:AZ, Cochise Co. 


Epargyreus clarus huachucal 


- Epargyreus clarus huachucal 15105Ell|Mexico:Sonora 
Epargyreus clarus profugus|11698|USA:NM,Grant Co. 

Epargyreus clarus profugusl 15038A11|HT|USA:NV,Clark Co. 


rtcbt:— Epargyreus clarus profugusl 15116H09|USA:AZ,Coconino Co 
Li o^ 24 Epargyreus clarus profugusl 15116H08|USA:AZ,Mohave Co. 

w*- Epargyreus clarus profugus] 15116H10|USA:AZ,Mohave Co. 

— Epargyreus clarus profugusl li636|USA:AZ,Green1ee Co. 


Epargyreus clarus profugus 11185 1 USA:NM,Otero Co. 


Epargyreus orizabal 14123C06|Belize:Cayo 
- Epargyreus cruza| 16107G12|Costa Rica:ACG|09-SRNP-73319 


Fig. 26. Epargyreus clams (blue) and huachuca (red) with outgroups. 


Nisoniades bromias (Godman & Salvin, 1894) is a species 
distinct from Nisoniades rubescens (Moschler, 1877) 

Considered a junior subjective synonym of 
Nisoniades rubescens (Moschler, 1877) (type 
locality Suriname), Pellicia bromias Godman & 

Salvin, 1894 (type locality Mexico, Guatemala, 

Costa Rica, Panama) reveals 1.2% (8 bp) 


— Nisoniades rubescens| 17116A07|M|Mexico:Tamaulipas| 1974 
Nisoniades rubescens|18059F12| Nicaragua 11967 

Nisoniades rubescens 118059F101 Colombial 1972 
- Nisoniades rubescens|4968|M|French Guianalold 

Nisoniades ephoral 18059G09|Panama| 1981 


Nisoniades mimas[l8059G12lPerui 1982 

Fig. 27. Nisoniades bromias (red) and N. rubescens (blue). 
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difference in COI barcode, which by itself is not large enough to draw definitive conclusions, but it 
prompted further investigation. The Fst/Gmin statistics computed on two pairs of specimens from distant 
localities (Fig. 27) were 0.27/0.05, suggesting that Nisoniades bromias (Godman & Salvin, 1894), 
reinstated status, is a distinct species. 


Pholisora crestar J. Scott & Davenport, 2017 is a species 
distinct from Pholisora catullus (Fabricius, 1793) 

Recently described as a subspecies on the basis of wing pattern differences (Scott et al. 2017), Pholisora 
catullus crestar J. Scott & Davenport, 2017 (type locality USA: California, Tulare Co.) was synonymized 
with Pholisora catullus (Fabricius, 1793) 

(type locality "Indiis", likely eastern US) 
by Pelham (2020). As a part of on-going 
genomic sequencing inventory of the 
primary type specimens of Hesperiidae, 
we obtained and analyzed whole genome shotgun reads of the holotype and two paratypes of crestar. 
Surprisingly, their comparison with P. catullus populations from several distant localities revealed 
prominent genetic differentiation (Fig. 28). Moreover, one of the crestar paratypes (from CA: Mono Co., 
NVG-17066H12, Fig. 28) apparently is P. catullus, not crestar. Fst/Gmin statistics for the crestar!catullus 
comparison are 0.34/0.014, suggesting distinctness of crestar as a species. Gene exchange between 
catullus and crestar (0.014) is very low (for conspecific populations it is typically above 0.1), strongly 
supporting reproductive isolation between these taxa. Peculiarities of COI barcode evolution in Pholisora 
have been reported previously by Pfeiler (2018) and COI barcodes of the crestar holotype and the catullus 
specimen from Texas (NVG-3990) differ by 1.7% (11 bp). Due to strong genetic differentiation, we 
suggest that Pholisora crestar J. Scott & Davenport, 2017, new status, is a species-level taxon. 


Pholisora catullus crestar|19083E02IHT|USA:CA,Tulare CO.I20131CSUC 
Pholisora catullus crestar|17066HUJPTjUSA:CA,Kern Co.|2010|CSUC 

Pholisora catullus catullus| 17066H12JPT of crestar|USA:CA.Mono Co.|2006|CSUC 
— Pholisora catullus catullus] 19042A05IMexico:Hidalgo|1964 
Pholisora catullus catu lusl4278[USA;IN.Montgomerv Co,|2015 
Pholisora catullus catullus|3990 USA:T/,Starr Co.f2015 

- Pholisora catullus catullusl i6107D02|USA:TXJeff Davis Co.|2016 

- Pholisora meiicanusl7976IUSA:NM,Colfax Co..Cimarron] 1989 

Fig. 28. Pholisora crestar (red) and P. catullus (blue). 


Carterocephalus mandan (W. H. Edwards, 1863) and 
Carterocephalus skada (W. H. Edwards, 1870) are species-level taxa and 
not subspecies of Carterocephalus palaemon (Pallas, 1771) 

Proposed as a species, Hesperia mandan W. H. Edwards, 1863 (type locality Canada: Manitoba) has 
mostly been considered a subspecies of Carterocephalus palaemon (Pallas, 1771) (type locality Russia: 
Samara Oblast) (Pelham 2008). Although it has been recently reinstated as a species (Pohl et al. 2010), 
this suggestion was not universally followed (Pelham 2020). Genomic comparison of the Old and New 
World palaemon-like populations reveals three clusters in the tree (Fig. 29): two corresponding to the 
abovementioned taxa (red and blue), and the third one (green) for Carterocephalus palaemon skada (W. 
H. Edwards, 1870) (type locality USA: AK, Kodiak) together with Carterocephalus palaemon magnus 
Mattoon & Tilden, 1998 (type locality USA: California, Sonoma Co.). The largest separation is observed 
between the Old World C. 
palaemon and the New 
World taxa. The Fst/Gmin 
statistics for the comparison 
of pairs of these 3 clusters 
(we call the green cluster by 
its oldest name: skada) are: 
palaemon vs. mandan 0.50/ 

0.016, palaemon vs. skada 
0.56/0.005, and mandan vs. skada 0.36/0.025. All of these numbers indicate strong genetic differentiation 
and very low gene exchange between clusters. Analysis of COI barcodes reveals an unusual situation. 
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— Carterocephalus palaemon mandan|15096D03lNT|Canada:Manitoba 
— Carterocephalus palaemon mandanl 16106E03|USA:MN.Lal.e of the Woods Co. 
Carterocephalus palaemon mandanl 16106E02|USA:MN,Lake of the Woods Co. 

o 62 - Carterocephalus palaemon skadal 17067C01|USA:AK,Fairbanks North Star 

—- Carterocephalus palaemon skada-cf (S. Rocky Mts. segregate)! 17067B11 |USA:WY,Teton Co. 


Carterocephalus palaemon magnus|PA069|USA:CA,Sierra C 
Carterocephalus palaemon magnus|PA0427|USA:CA,Mendocino Co. 

- Carterocephalus palaemon magnus|17109A05|HT|USA:CA,Sonoma Co. 
Carterocephalus palaemon magnus| 17067C03|USA:CA,Sonoma Co. 

Carterocephalus palaemon magnus|17067C02|USA:CA,Sonoma Co. 

- Carterocephalus palaemon palaemon (=ab. carrueli)| 18079B07|HT|France 


Carterocephalus palaemon palaemonl 17067B08|Germany:Bayern 
Carterocephalus palaemon palaemonl 18092F02|Poland 
_srocephalus palaemon palaemonl 17067B09|Russia:Orenburg Prov. 

Carterocephalus silvicola[18038G04JRussia:Buryatia 
— Carterocephalus*. 


Carterocephalus silvicola118123All|Sweden 


Fig. 29. Carterocephalus mandan (red), skada (green) and palaemon (blue). 



































First, barcode difference of palaemon vs. mandan is the same as palaemon vs. skada : 1.5% (10 bp). 
Second, mandan barcodes are not much different from skada (0.3%, 2 bp), which can be explained by 
introgression, but the mandan neotype and magnus holotype exhibit larger difference of 0.76% (5 bp) 
between them. Third, for comparison, two widely sympatric Old World species C. palaemon and 
Carterocephalus sdvicola (Meigen, 1829) exhibit Fst/Gmin of 0.72/0.0006 (indicating very strong 
isolation), but barcode difference between them is only 0.6% (4 bp). A number of similar instances of 
distinct butterfly species not strongly different in their barcodes have been documented (Burns et al. 2008; 
Cong et al. 2017), thus barcode differences and similarities cannot be considered separately from all other 
evidence. In summary, we suggest to reinstate Carterocephalus mandan (W. H. Edwards, 1863) and 
Carterocephalus skada (W. H. Edwards, 1870) as species, and additionally propose the following revised 
combinations: Carterocephalus mandan mesapano (Scudder, 1868) and Carterocephalus skada magnus 
Mattoon & Tilden, 1998. 


Amblyscirtes arizonae H. Freeman, 1993 is a species 
distinct from Amblyscirtes elissa Godman, 1900 

Described as a subspecies of Amblyscirtes elissa Godman, 1900 (type locality Mexico: Guerrero), A. e. 
arizonae H. Freeman, 1993 (type locality USA: Arizona, Santa Cruz Co.) has not caused much attention 
being invariably kept as a subspecies 
(Pelham 2008). We sequenced primary 
type specimens of both A. e. arizonae 
and A. elissa and compared them with 
specimens from other localities. The 
genomic tree revealed two prominent 
clusters with separation similar to that 
between A. belli H. Freeman, 1941 and 
A. celia Skinner, 1895, and between A. 

Carolina (Skinner, 1892) and A. reversa 
F. Jones, 1926 (Fig. 30). The Fst/Gmin 
statistics for them were 0.60/0.002, implying strong genetic differentiation and virtually no gene exchange 
between these taxa. The COI barcodes of A. elissa and A. e. arizonae primary type specimens differ by 
2.6% (17 bp). For these reasons, we suggest that Amblyscirtes arizonae H. Freeman, 1993 is a distinct 
species, new status. 

Megathymus violae D. Stallings & Turner, 1956 is a species 
distinct from Megathymus ursus Poling, 1902 

Initially proposed as a species, Megathymus violae D. Stallings & Turner, 1956 (type locality USA: New 
Mexico, Eddy Co.) was placed as a 
subspecies of Megathymus ursus 
Poling, 1902 (type locality USA: 

Arizona, Pinal Co.) by dos Passos 
(1960). Genomic comparison of 
specimens of both taxa across their 
ranges, including the holotypes, 
revealed their prominent separation 
in the tree with the distance close to 
that between Megathymus yuccae 
(Boisduval & Le Conte, [1837]) and 




Amblyscirtes elissa arizonael 15096A05|HT|M|USA:AZ,Santa Cruz Co. 

• Amblyscirtes elissa arizonael 19064B12|M|USA:AZ,Santa Cruz Co.11986 
Amblyscirtes elissa arizonae|9685|F|USA:AZ,Santa Cruz Co.|2017 
Amblyscirtes elissa arizonael 18063F06|M|Mexico:Sonora| 1988 

- Amblyscirtes elissa elissa118083E05|ST|M|Mexico:Guerrero 
Amblyscirtes elissa elissal 18063F12|M Mexico:Chiapas| 1989 

- Amblyscirtes elissa elissal 19042G02|M|Mexico:Guerrero| 1951 

- Amblyscirtes elissa elissal 19042G03|M|Mexico:Morelos| 1942 
..Amblyscirtes celia|7644|M|Mexico:Tamaulipas / Ciudad Mantel 1973 
Amblyscirtes celia[4501|M|USA:TX,Cameron Co,E of Brownsviile|2015 

0 , T Amblyscirtes celia| 15096H02|ST|M|USA:TX,Comal Co.|CMNH 
Amblyscirtes celia| 1671|USA:TX ; Travis Co.,Platt Lane|2011 
Amblyscirtes belli| 18025E02|HT|M|USA:TX,Dallas Co. 
f-r Amblyscirtes belli|3699|M|USA:TX,Dallas Co,Dallas|2015 
L Amblyscirtes belli| 1240|USA:TX,Brazos Co.,Lick Creek Park|2013 

i—- Amblyscirtes cardinal 15095G08|ST|USA:VA,Richmond Co. 

Lr 2 *?— Amblyscirtes cardinal 10672|M|USA:VA,Chesapeake| 1970 

“- Amblyscirtes carolina| 18012C05|M|USA:NC,Richmond Co.11984 

Amblyscirtes reversal 15096C04|HT|USA:VA,Nansemond Co. 

0 T. Amblyscirtes reversal 18012C03|F|USA:NC,Brunswick Co.11995 
- Amblyscirtes reversal 18012C06|M|USA:NC,Cumberland Co.11984 
Amblyscirtes reversal 10673|M|USA:VA,Virqinia Beach|1971 


I - ZJ 


Fig. 30. Amblyscirtes arizonae (red), A. elissa (blue) and others. 
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Megathymus beulahae D. Stallings & J. Turner, 1958 (Fig. 31). The Fst/Gmin statistics for comparison of 
ursus and violae groups are 0.56/0.001 (note close to 0 gene exchange between these taxa). The COI 
barcodes of the M. ursus and M. violae holotypes differ by 1.8% (12 bp). For these reasons, we reinstate 
Megathymus violae D. Stallings & Turner, 1956 as a species-level taxon. 


Discussion: genomic trees, branch lengths and genera 

Near the end, coming back to the Introduction, we elaborate on and illustrate the reasons behind the 
classification decisions that we have chosen to make about genera. Traditionally, species were grouped 
into genera by phenotypic characters. For butterflies, these were mostly wing patterns and shapes, and 
genitalic morphology. When differences in these phenotypic aspects were deemed to be significant 
enough according to a subjective opinion of an individual researcher, they formed a basis for defining a 
genus. This system served its purpose until a consensus opinion was formed among taxonomists that each 
genus should be monophyletic. It is exceedingly difficult to predict monophyletic taxa from their 
phenotypes, and DNA-based phylogenetic trees provide the most reliable inference of monophyletic 
groups. Therefore, genera should be defined using phylogenetic trees constructed from DNA sequences. 

Each individual feature of an organism can experience rapid evolution and fool researchers into 
making incorrect classification decisions. Genitalia that are commonly used in Lepidoptera classification 
are prone to such rapid changes as well. For instance, Steinhauser (1989) proposed a genus Thessia on the 
basis of unique shape of genitalic valvae. However, even a very short, 654 base pair region of DNA, such 
as the COI barcode, reveals the paraphyly of Achalarus Scudder, 1872 (as it was circumscribed at that 
time) with respect to Thessia (Pfeiler et al. 2016), suggesting that the unique valva is a result of 
accelerated evolution within Achalarus rather than a character originated after Thessia and Achalarus 
have (supposedly) diverged from each other. Therefore, a decision to erect the genus Thessia was a 
mistake, because Thessia is a subclade within (as it was then defined) Achalarus. Nevertheless, the 
barcode DNA region itself is a single feature, and as any other such feature, can experience evolutionary 
irregularities. To reduce such mistakes, it is better to use information from as many features as feasible. 
Complete genomes offer the ultimate DNA dataset for classification decisions. Genomic analysis suggests 
that Achalarus itself is a junior subjective synonym of the subgenus Thoiybes Scudder, 1872, and Thessia 
is actually a junior subjective synonym of the subgenus Murgaria E. Watson, 1893 (Li et al. 2019). 

Genomic trees summarize integral information about the entire organism, not just some of its 
features. For this reason, we use them to make decisions about classification of genera. Here, we explain 
how we arrive to these decisions using examples from this work and our previous publication (Zhang et 
al. 2019c). A maximum likelihood tree constructed using IQ-TREE program (model GTR+I+G) (Minh et 
al. 2020) from concatenated protein-coding regions of nuclear genomes is shown in Fig. 32. To best 
follow our logic, a reader may close the tree on the right (Fig. 32b, the final result) and look only at the 
tree on the left (Fig. 32a), which is the same as the tree on the right, but without the final results being 
marked in order not to bias the reader. This tree was constructed without assuming a molecular clock and 
reveals differences in evolutionary rates between species: i.e., species names are placed at difference 
distance from the left side of the page (=from the root of the tree). We see that Emesis evolved the fastest 
(the farthest from the left), and Ephyriades Hiibner, [1819] evolved the slowest (closest to the left). In a 
tree, only horizontal (left-to-right) distances matter. Vertical (top to bottom) distances are arbitrary and 
are set to place species names evenly along vertical dimension, so that the names do not overlap and are 
not too far away from each other to save space. 

Tree branches have different lengths. Again, only horizontal branches have evolutionary meaning, 
and vertical lines in the tree are set to avoid overlap of names and to connect branches to nodes. The 
length of a horizontal branch is proportional to the number of estimated changes in DNA (=fixed 
mutations) that happened along the branch. The tree has a scale bar near the bottom (Fig. 32). The length 
of that bar, as indicated, corresponds to 6 changes per 100 base pairs (=0.06, or 6%). Using this bar, we 
can measure evolutionary distances between taxa in DNA changes. Long branches correspond to many 
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Fig. 32. Delineating genera using a tree constructed from nuclear genomes: unmarked (a) and colored by suggested genera (b). 


changes in genomic DNA. Short branches correspond to few changes in genomic DNA. Because larger 
number of DNA changes are expected to result in larger number of phenotypic changes, longer branches 
correspond to more phenotypic changes on average. These are integral changes and some of them may be 
in genitalia, others may be in caterpillar morphology. Regardless of where these changes are, longer 
branches are more important than shorter branches. In addition to larger number of changes, longer 
branches are also more reliable and support clades that are more likely to be correct. The statistical 
reliability of every clade is indicated by a number next to each node. This number is a fraction of trees 
(out of 100 trees constructed from various subsets of genomic segments) that contain this node, e.g. a 
genome was divided into 100 segments and each segment was used to generate a tree. If a particular node 
is present in all 100 trees, the number by that node is 1. Therefore, this number measures consistency 
between trees constructed from different partitions of the data. If every DNA segment supports a clade, it 
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has a number 1 next to it. If 94 out of 100 segments support the clade, the number is 0.94. 

A genus should be a prominent, major clade in the tree that is above species level and below tribe 
and sub tribe levels. Phenotypic features are difficult to quantify, and due to the possibly uneven speed of 
evolution, it is a challenge to detennine which phenotypic changes correspond to major clades. Total 
genomic changes can be used as a yardstick to quantify each clade. The number of total genomic changes 
is proportional to branch lengths in genomic trees (Fig. 32a). Therefore, the task of identifying genera 
may be viewed as a task of identifying prominent (i.e. supported by longer branches compared to 
surrounding branches) clades in genomic trees that on average correspond to how genera are defined 
currently (to avoid unnecessary taxonomic changes). Additionally, we believe that each genus should not 
be very different from another genus in terms of genetic differentiation of species placed in a genus, i.e. 
genera could be defined consistently, so that genera correspond to clades of approximately the same 
differentiation within. Defined consistently, the genus becomes a level (as meant by this word) of a 
classification instead of several varying levels, i.e., we can expect a genus to be a group of species bearing 
about the same relatedness among them as that in other genera. It would seem unnatural if one 
phylogenetic group is oversplit into genera, i.e. genera in that group correspond to very closely related 
species, but another group is undersplit, and genera in it correspond to species that are only distantly 
related. The measure of closeness as we use it, is overall genomic divergence. 

Looking at the clade of Hesperiidae at the top of the tree (Fig. 32a) we see three major clades, not 
two and not four. The first clade is Ephyriades and is sister to all other taxa. Then all others split into two 
clades of similar genetic differentiation within each clade. We see that each of these clades resembles a 
tight bush or a comb, rather than an evenly bifurcating tree, i.e. the internal branches in either clade are 
much shorter than a branch that supports the entire clade. The clade with Gesta bifurcates into two 
subclades, one consists of Gesta sensu stricto (s. s.). Species from the other subclade were called 
"Erynnis" previously (and are called Erynnis in the tree to facilitate communication): it is a subgenus 
Erynnides Burns, 1964 (type species Nisoniades propertius Scudder & Burgess, 1870). If we consider 
these two subclades to be major clades, then the Hesperiidae tree would consist of four major clades 
(. Ephyriades, Erynnis s. s., Erynnides and Gesta). However, the branches supporting the two subclades 
(Erynnides and Gesta) are nearly three times shorter than the branches supporting the clades Erynnis s. s. 
and a clade combining Erynnides with Gesta. Therefore, the Hesperiidae subtree should not be partitioned 
into four major clades, because two of these clades ( Erynnides , Gesta) would be minor compared to the 
other two, and more importantly, compared to the clade combining Erynnides with Gesta. 

The remaining alternative to a three-clade partitioning would be a two major clade partition, where 
Erynnis s. s., Erynnides and Gesta are all joined together into Erynnis sensu lato (s. 1.) The branch 
supporting this clade is only slightly shorter than the branch supporting Erynnis s. s., and therefore this 
clade is rather prominent in the tree. We reject this solution for the two reasons. First, Erynnis s. 1. is not a 
homogenous group of species, which we think a genus should be, i.e. the Erynnis s. 1. clade does not look 
like a bush or a comb. Instead, it splits into two major clades: Erynnis s. s. and Erynnides + Gesta, (we 
call this clade Gesta s. 1. from now on) each of which individually looks more like a comb than when they 
are combined. In other words, Erynnis s. 1. itself is composed of two major clades, and does not represent 
a single group of species, but two major groups of species. 

The second reason stems from consistency between different genera, i.e. an idea that different genera 
should represent the same level in the classification (Fig. 33). Being a level, genera should be groups of 
species with comparable divergence within each genus. In this tree (Fig. 32a), where all branches are to 
scale, we can compare divergence between Erynnis s. s. and Gesta s. 1. to the divergence in Nymphalidae 
previously placed in genera Aglais, Polygonia, Nymphalis, and Vanessa. These two subtrees (Erynnis and 
Vanessa) are illustrated in Fig. 33. Genetic differentiation of a clade is proportional to the average 
distance (average sum of branch lengths) from the last common ancestor of the clade (=node that supports 
the entire clade) to the leaves (=species) in the clade. In other words, it is a linear distance (in horizontal 
dimension) from the base of the clade to the tips of the tree. On the one hand, we see that Polygonia 
divergence is rather small, perhaps comparable to the divergence of the Erynnides subclade with horatius 
and Juvenalis, and definitely smaller than the divergence within either Erynnis s. s., or Gesta s. 1. On the 
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other hand, the divergence of Erynnis s. 1. is larger than 
the divergence of Aglais, Polygonia, Nymphalis and 
Vanessa combined. Therefore, having Erynnis s. 1. as a 
genus is inconsistent with having Polygonia as a genus: 
these two groups represent different levels in the 
classification. Coming back to Nymphalidae, we see that 
branches supporting Aglais, Polygonia and Nymphalis 
individually are much shorter than the branches 
supporting Erynnis s. s. or Gesta s. 1. Only the branch 
supporting Vanessa is somewhat comparable, although 
shorter. However, the branch supporting the first three 
clades together {Nymphalis s. 1.) is more prominent and 
is about the same as the branch supporting Vanessa. 

In summary, Eiynnis s. 1. is comparable to Vanessa 
s. 1. {Nymphalis s. 1. + Vanessa s. s.). A system of two 
genera {Erynnis s. s. and Gesta s. 1.) is comparable to two 
genera Nymphalis s. 1. and Vanessa s. s. We attempt to 
choose an internally consistent solution that agrees the 
most with how these species are assigned to genera in the 
current classification. Therefore, we choose the 2-genus 
solution for both of these cases, as shown in Figs. 32b 
(colored clades E: Erynnis, G: Gesta, N: Nymphalis and 
V: Vanessa) and 33 (shaded clades). These four genera 
represent a similar level in the classification and correlate 
with the current classification of these butterflies. The 
choice of Erynnis s. 1. would correspond to a consistent 
choice of joining all four Nymphalidae genera in 
Vanessa, which may represent too much of a lump and 
more name changes (Fig. 33). 

Another point is that genetic differentiation can be 
used to estimate divergence times of these clades through 
the tree rescaling and calibration with fossils (primary 
calibration) (Chazot et al. 2019) or other time-calibrated trees (secondary calibration) (Zhang et al. 
2019a). As we have seen in Hesperiidae (Li et al. 2019), the genus level typically corresponds to 
divergence between 10 and 15 million years ago (Mya). Divergence of Erynnis s. 1. was estimated to be 
about 27 Mya, which is larger than the divergence between Vanessa s. s. and Nymphalis s. 1., at about 22 
Mya (Zhang et al. 2019d). However, divergences within Gesta s. 1. (~ 16 Mya), Vanessa s. s. (~ 16 Mya) 
and Nymphalis s. 1. (-14 Mya) (Zhang et al. 2019d) are very much comparable to each other, and these 
genera represent groups of about the same level. It should be noted that the divergence times are only 
approximate, should be considered with caution, and may have errors of possibly up to 50%, especially in 
groups with large differences in evolutionary rates. However, the relative comparison of divergence times 
estimated within the same tree using the same method is expected to be more accurate. 

Finally, a question arises about how these considerations of trees, branch lengths, divergence and 
geological times correlate with genera definition based on phenotypic characters. Because phenotypic 
characters are encoded by the genotype, longer branches in the tree that correspond to more changes in a 
genotype (these are integral genomic trees, not based on several gene markers) should translate to more 
changes in the phenotype. We advocate a method to delineate genera from genomic trees first, and then 
come back to phenotypic analysis to find the phenotypic characters that correspond to these genera. In the 
case of Erynnis and Gesta, the retrospective inspection of morphological characters yields substantial 
differences in male genitalia that have been noted and illustrated previously (Evans 1953; Bums 1964). 
The uncus is asymmetric, terminally broad in Gesta, but is symmetric, extending into a "beak" in Erynnis. 
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The valvae are strongly asymmetric with at least one extended harpe in Gesta, but are more symmetric 
with shorter harpes in Erynnis. Other differences are stated in the diagnosis of Erynnides by Burns (1964). 

Comparing the clades of other groups in Fig. 32a with Erynnis/Gesta and Nymphalis/ Vanessa we see 
that divergence within Speyeria and Roeberella (a clade containing R. clavus and with Apodemia 
hypoglauca at its base), and divergence between Hypaurotis, Favonius and Elabrodais is much smaller 
than that in the groups we define as genera. We also see that the colored clades (with letters denoting 
corresponding genera by each clade) in Fig. 32b are more or less equivalent to each other in terms of 
genetic differentiation (distance from the base of the clade to its tips) and prominence (length of the 
branch supporting the clade). For these reasons, we suggest that these clades can be treated as genera: 
they are prominent, consistent, and reasonably well correspond to how genera have been defined 
previously. The changes we suggest combine some more compact in terms of genetic (and phenotypic) 
differentiation genera into more internally diverse genera that become more consistent with the 
differentiation within many classic genera such as Emesis, Ministrymon, Vanessa, and Boloria. 
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